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SECTION  M 
INTRODUCTION 


This  report  presents  final  results  from  a  FY -89  VLF/LF  Data  Analysis  task  initiated 
in  support  of  the  VLF/LF  propagation  modeling  program  being  developed  by  the  Naval 
Ocean  Systems  Center  (NOSC).  Specifically,  this  paper  gives  the  results  of  comparisons 
between  VLF  signal  measurements  recorded  during  nighttime  hours  and  NOSC  model 
predictions  for  similar  nighttime  propagation  paths.  Measured  data  was  recorded  in 
1985  and  1986  aboard  the  merchant  ship  GTS  Callaghan,  and  daytime  propagation  com¬ 
parisons  against  that  data  were  documented  in  the  NOSC  report  entitled  “Evaluation 
of  the  NAVOCEANSYSCEN  Long  Wavelength  Propagation  Capability  Using  VLF  Data 
Collected  Aboard  a  Ship”,  01  November  1987  (termed  herein  the  Callaghan  Report). 
The  present  report  is  offered  in  satisfaction  of  CDRL  A002,  under  D.  0.  0046,  Contract 
N66001-88-D-0033. 

For  many  years,  NOSC  has  developed  computer  programs  and  models  for  the  analysis 
of  VLF/LF  propagation,  and  has  applied  this  capability  to  a  number  of  applications  for 
the  Navy  and  for  other  DoD  agencies.  These  past  efforts  have  utilized  separate  pro¬ 
grams  for  various  parts  of  the  calculations  and  have  generally  been  performed  by 
experienced  NOSC  personnel  capable  of  judging  the  correctness  of  resulting  outputs. 

More  recently,  NOSC  has  developed  a  new  approach  to  longwave  computations, 
whereby  many  parameter  selections  are  automated,  permitting  reliable  use  by  less  skilled 
operators.  This  integrated  program  approach,  called  the  LongWave  Propagation  Capa¬ 
bility  (LWPC),  is  now  in  the  initial  release  process  and  is  expected  to  find  broad 
application  within  the  defense  community. 

Verification  of  the  program’s  integrity  through  a  continuing  process  of  comparison 
with  existing  and  future  VLF/LF  field  measurements  is  needed.  Such  comparative  anal¬ 
yses  guide  the  development  of  the  ionospheric  model  used  in  the  LWPC,  which  must 
be  fine-tuned  to  assure  maximum  accuracy  when  applied  to  increasingly  complex  world¬ 
wide  scenarios. 

In  an  earlier  effort  by  NOSC,  VLF/LF  data  recorded  on-board  the  merchant  ship 
GTS  Callaghan  was  reduced  and  aggregated  using  an  HP-9000  computer,  and  was  sub¬ 
sequently  compared  with  predicted  data  using  the  LWPC  model.  The  comparison  was 
limited  to  four  VLF  frequencies,  and  to  davtime-only  conditions.  In  the  present  task, 
similar  analyses  of  nighttime  data  are  conducted,  again  limited  to  the  VLF  band. 


1 


In  the  sections  that  follow,  this  report  summarizes  certain  adjustments  recommended 
for  the  LWPC  ionospheric  model  as  the  result  of  the  comparisons,  and  gives  some  con¬ 
clusions  about  analysis  techniques  that  could  be  used  in  future  comparisons.  The  main 
body  of  the  report  is  given  in  section  II,  where  the  objectives  and  methodology  of  the 
analysis  are  stated,  along  with  discussions  of  the  data  reduction  and  propagation  pre¬ 
diction  computer  programs  utilized  in  the  analysis.  Analysis  of  the  Callaghan  nighttime 
data  is  discussed  in  section  11-5,  and  sample  comparisons  of  predicted  and  measured  data 
are  presented  in  section  II-6.  Statistical  evaluations  of  the  data  are  shown  in  section 
II-7  along  with  summary  charts  that  provide  a  technique  for  optimized  selection  of  ion¬ 
ospheric  profiles.  Final  conclusions  and  recommendations  are  discussed  in  section  II-8. 
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SECTION  1-2 

SUMMARY  AND  CONCLUSIONS 


This  paper  documents  a  VLF/LF  Data  Analysis  task  by  the  Naval  Ocean  Systems 
Center  to  improve  the  modeling  of  the  nighttime  ionosphere  when  making  propagation 
predictions  with  the  LWPC  computer  program.  The  task  utilizes  an  extensive  database 
of  VLF  measured  data  recorded  during  the  1985-1986  trips  of  the  merchant  ship  GTS 
Callaghan  in  the  North  Atlantic  area.  By  constraining  the  Callaghan  data  to  those  pe¬ 
riods  when  both  the  ship  and  the  distant  transmitters  were  in  time  zones  consistent  with 
all-nighttime  propagation,  and  by  eliminating  data  from  trips  outside  the  principal  area 
of  inteiest,  an  aggregated  set  of  recorded  data  was  assembled  for  each  frequency  of  con¬ 
cern.  Four  frequencies  were  examined:  16.0,  19.0,  21.4  and  24.0  kHz. 

Recorded  data  sets  were  graphed  as  signal  vs.  distance  plots,  computing  distance  from 
the  transmitter  for  each  ship’s  location.  The  LWPC  program  was  then  utilized  to  com¬ 
pute  signal  vs.  distance  along  a  typical  path  in  the  same  ocean  area,  and  the  predicted 
and  recorded  data  were  compared.  By  changing  the  LWPC  parameters/!  and  h',  which 
are  exponential  relationships  describing  ionospheric  electron  density  and  collision  fre¬ 
quency  profiles  that  vary  with  height,  different  propagation  predictions  were  compared 
with  the  recorded  data  until  a  best  fit  was  obtained. 

In  all,  some  265  separate  LWPC  runs  were  compared  with  measured  data  at  four 
frequencies,  using  manual  overlays  on  a  light  table  to  find  the  best  match  between  pre¬ 
dicted  and  measured  curves.  Although  this  technique  enabled  the  selection  of  improved 
/3-h'  values,  the  method  was  labor  intensive  and  required  subjective  judgments  that  be¬ 
came  difficult  at  the  higher  frequencies.  However,  best-fit  ^-h'  values  for  each 
frequency  were  selected  by  that  method,  and  a  final  confirmation  was  made  using  a  sta¬ 
tistical  evaluation  technique  wherein  average  values  of  the  means  and  standard 
deviations  of  the  difference  between  recorded  and  predicted  signal  levels  were  graphed. 

The  present  LWPC  default  values  for  nighttime  /3  and  h'  are  0.30-0.80  km-1  and  87 
km,  respectively,  (fi  is  entered  as  a  range  corresponding  to  the  frequency  range  of  10-60 
kHz;  LWPC  then  finds  (i  for  a  specific  frequency  by  interpolation.)  The  analysis  effort 
determined  that  a  0  range  of  0.40-0.90,  and  an  h'  of  85  provided  a  better  fit  between 
recorded  and  predicted  data,  with  the  improvement  becoming  quite  significant  at  the 
lower  frequencies.  A  preliminary  recommendation  for  LWPC  /J-h'  profiles  is  shown 
in  figure  1. 
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For  the  higher  frequencies  of  21.4  and  24.0  kHz,  propagation  paths  from  the  CONUS 
sources  of  these  signals  fell  above  the  70°  dip  angle  line  for  a  portion  of  the  paths  (see 
figure  3).  Since  the  LWPC  program  adjusts  the  default  p-h'  profiles  for  polar  cap 
(PCAP)  effects  on  paths  traversing  areas  above  the  70°  line,  and  since  part  of  the  task 
was  to  determine  the  validity  of  the  PCAP  default  as  well  as  that  of  the  /9-h'  profiles, 
the  analysis  of  the  higher  frequencies  became  more  complex.  To  examine  PCAP  effects, 
duplicate  LWPC  runs  were  made.  Normal  runs  utilized  the  default  PCAP  setting  of  70- 
74,  while  excursion  runs  set  the  PCAP  range  at  80-84,  which  had  the  effect  of  moving 
the  dip  angle  line  well  above  the  propagation  paths  of  interest,  thus  preventing  automatic 
LWPC  profile  modification.  The  results  of  these  tests  showed  better  data  fits  for  the 
excursion  settings,  indicating  that  the  default  PCAP  values  may  require  adjustment  at 
the  higher  VLF  frequencies.  Results  for  both  settings  are  shown  in  figure  1. 

It  was  concluded  that  the  present  LWPC  default  settings  for  nighttime  ionospheric 
profiles  should  be  adjusted.  Preliminary  profile  recommendations  are  to  increase  the 
P  range  by  0.10  km*1  and  to  decrease  h'  by  2  km.  The  default  PCAP  settings  of  70-74 
may  also  require  adjustment,  but  further  data  analysis  is  required  before  further  sug¬ 
gestions  may  be  made  in  this  area.  The  following  recommendations  are  made: 

•  LWPC  profile  optimization  should  be  mechanized  by  modifying  the  statis¬ 
tical  graphing  program  to  automatically  determine  the  average  of  mean  and 
standard  deviation  values  of  recorded-minus-predicted  data,  for  a  number 
of  distance  ranges.  Resulting  averages  should  be  plotted  as  a  function  of 
p  and  h'  as  a  means  of  finding  truly  optimized  profiles  for  each  frequency. 

•  The  optimization  of  LWPC  nighttime  profiles  should  be  extended  over  a 
much  wider  range  of  p-h'.  This  larger  effort  would  be  made  feasible  by 
the  automatic  statistical  graphing  technique  suggested  above. 

•  The  more  extensive  LWPC  runs  should  be  facilitated  by  utilizing  the  newer, 
machine-independent  LWPC,  Version  1.0,  which  will  permit  faster  runs  and 
eliminate  the  need  for  VAX- VMS  procedures. 

•  Daytime  profile  optimization  of  the  LWPC  as  determined  in  the  Callaghan 
Report  should  be  re-examined,  using  automated  statistical  profile  compar¬ 
isons  over  a  wider  p- h'  range. 

•  Both  daytime  and  nighttime  comparisons  should  be  extended  to  higher  fre¬ 
quencies  as  measured  data  becomes  available.  At  a  minimum,  the  LF 
channels  of  the  Callaghan  database  should  be  examined. 
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FREQUENCY  (kHz) 


FREQUENCY  (kHz) 

-H-  Preliminary  Q  With  —  LWPC  Default 

Recommendation  PCAP=80,84  Profile 

With  PCAP=70,74 


Figure  1.  Preliminary  recommendation, 

revised  /3— h*  ranges  for  LWPC  nighttime  propagation 
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SECTION  IM 
OBJECTIVES 

The  overall  task  objective  is  to  provide  an  in-depth  comparative  analysis  between  ex¬ 
isting  NOSC  VLF/LF  field-strength  measurements  and  signal-strength  predictions  as 
generated  by  the  LWPC  computer  program.  A  principal  goal  is  to  develop  recommen¬ 
dations  for  the  specification  of  refined  LWPC  ionospheric  (/9-h')  profiles,  so  that  future 
LWPC  predictions  agree  more  closely  with  observed  data.  In  this  task,  the  analysis  is 
constrained  to  VLF-band  frequencies  and  nighttime  propagation  only. 

An  interim  report*  documented  the  following  task  elements: 


•  Review  documentation  and  computer  programs  which  comprise  the 
LWPC.  Perform  LWPC  test  computations  on  contractor  computers 
and  compare  results  with  similar  calculations  supplied  by  NOSC. 


•  Review  NOSC-supplied  LWPC  predictions  as  given  in  the  1987  Cal¬ 
laghan  report,  and  compare  with  daytime  propagation  data  from  the 
GTS  Callaghan  measurements.  Develop  a  standard  methodology 
for  graphic  comparison  of  Callaghan  data  with  new  LWPC  calcula¬ 
tions  generated  within  the  task.  Compare  with  1987  Callaghan 
report  results  in  order  to  verify  analysis  methodology. 

The  present  report  documents  results  from  the  principal  task  elements,  which  are: 


•  Perform  LWPC  calculations  for  nighttime  cases  and  compare  with 
measurements  furnished  by  NOSC. 

•  Provide  comparative  analyses  of  measured  data  and  LWPC  predic¬ 
tions.  Provide  recommendations  for  refinements  in  LWPC 
ionospheric  model  to  better  align  measured  and  computed  data. 


•"Comparative  Analysis  -  Callaghan  Data”,  TCS  Report  TP-89-177,  31  May  1989. 
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SECTION  II-2 
METHODOLOGY 

The  methodology  of  the  VLF/LF  Data  Analysis  task  was  constrained  by  that  of  the 
Callaghan  Report;  every  effort  is  made  to  utilize  common  methods  and  data  reduction 
techniques  to  assure  a  continuity  of  results.  The  basic  technique  consists  of  comparing 
plotted  data  from  the  Callaghan  recorded  measurements  with  similar  plots  of  LWPC  sig¬ 
nal  predictions.  The  LWPC  ionospheric  profile  parameters  p  and  h'  are  then  adjusted 
as  necessary  to  obtain  a  best  fit  between  the  two  curves. 

Data-reduction  programs  supplied  by  NOSC  were  modified  and  used  to  produce 
curves  of  three  primary  types: 


•  Amplitude  vs.  distance  plots  of  raw  data  from  each  trip  of  the  GTS 
Callaghan. 

•  Amplitude  vs.  distance  plots  of  aggregated  data  from  trips  grouped 
by  number,  season  and/or  time-of-day. 

•  Statistical  analysis  charts  which  permit  quantitative  evaluations  of 
recorded  vs.  predicted  data  comparisons. 

Data  comparisons  and  best-fit  adjustments  of  P  and  h'  were  made  by  overlaying  pre¬ 
dicted-data  plots  with  those  of  aggregated  recorded  data,  using  the  analyst’s  best 
judgment.  Statistical  analysis  charts  were  then  produced  to  quantify  the  best-fit  choices. 

Analysis  of  a  typical  new  case  (a  new  transmitter  and/or  frequency)  followed  these 
procedures: 

1.  Choose  a  transmitter,  frequency,  time-of-day  and  season. 

Produce  plots  of  aggregated  recorded  data,  selecting  a 
group  of  trips  by  number  from  summary  lists  of  ship’s 
movements.  Figure  2  illustrates  a  group  of  Callaghan 
trips,  plotted  on  a  gnomonic  projection  map. 

Individual-trip  data  plots  may  be  consulted  as  necessary 
in  making  the  selection.  Alternately,  data  from  all  50 
trips  may  be  used  in  the  aggregated  plots. 
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Figure  2.  Map  showing  Callaghan  trips  selected  for  21.4  &  24.0  kHz, 
with  70°  dip  angle  limit  shown  as  wide  line. 


2.  Choose  a  propagation  path  from  the  desired  transmitter 
that  best  approximates  the  North  Atlantic  ocean  area  tra¬ 
versed  by  the  selected  ship  movements.  Figure  3* 
illustrates  the  four  propagation  paths  utilized  in  the  pres¬ 
ent  task  effort.  The  wide  line  denotes  the  position  of  the 
70°  dip  angle  limit. 

3.  Using  the  LWPC  program,  compute  signal  vs.  distance 
data  for  the  selected  path,  transmit  frequency  and  radi¬ 
ated  power  level.  Produce  a  plot  for  comparison  with  the 
aggregated  recorded  data  plot. 

4.  Iterate  the  LWPC  calculation  as  necessary  over  a  range 
of  ionospheric  profiles.  Select  best-fit  comparisons  with 
the  recorded  data  plot.  It  may  be  necessary  to  select  sev¬ 
eral  predicted  data  plots  for  later  selection  during  the 
statistical  analysis  process. 

5.  Produce  a  statistical  analysis  of  each  best-fit  selection  as 
a  plot  of  signal  histograms  vs.  distance.  Compute  the  av¬ 
erage  of  the  mean  and  standard  deviation  values  for  each 
applicable  distance.  Plot  the  averages  as  a  function  of  0 
and/or  h'  for  final  best-fit  selection. 


Because  certain  of  the  Callaghan  data  trip  files  were  found  to  be  unsuitable,  it  was 
necessary  to  make  a  preliminary  selection  of  recorded  data  to  be  aggregated.  Section 
II-5  discusses  the  factors  considered  in  this  selection. 


•The  map  of  figure  3  was  made  using  four  stacked  runs  of  the  NOSC  program  WMAP.  TCS 
has  modified  WMAP  to  permit  spooling  to  a  file,  with  later  off-line  printing  on  a  laser  printer. 
The  map  of  figure  2  was  generated  in  1  min.  30  sec.  on  the  AST-386  computer,  and  was  printed 
in  1  min.  14  sec. 
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Figure  3.  Propagation  paths  used  for  signal-strength  calculation 
as  a  function  of  distance,  with  70°  dip  angle  line  show 


SECTION  IIO 

DATA  REDUCTION  COMPUTER  PROGRAMS 

The  analysis  project  utilized  a  number  of  data  reduction  programs  and  supporting 
files  that  were  employed  by  NOSC  in  the  previous  Callaghan  data  analysis.  In  all,  some 
80  files  were  received  in  formats  suitable  for  use  on  Hewlett-Packard  Model  9000  com¬ 
puters.  The  project  examined  all  files  and  converted  a  number  of  the  Fortran  programs 
for  execution  on  MS-DOS  computers. 

The  following  NOSC  codes  were  converted  and  tested  for  routine  use  by  the  project: 

•  FORMIT  Forms  a  binary  data  file  for  a  specific  trip  and  frequency.  Gen¬ 

erates  an  optional  plot  of  signal  vs.  distance. 

•  PRHEAD  Reads  first  and  last  records  of  FORMIT-produced  binary  files. 

•  PLOT.DIS  Produces  aggregated  signal  data  at  a  single  frequency  of  selec¬ 

ted  trips  for  daytime,  nighttime  or  day-night  transition  periods. 
Plots  resulting  signal  data  as  a  function  of  distance.  Requires 
individual  trip  files  from  FORMIT  as  inputs. 

•  TRIPSTAT  Produces  a  statistical  analysis  of  the  difference  between  aggre¬ 

gated  recorded  data  and  predicted  data  vs.  distance.  Samples 
at  user-selected  intervals  (now  125  km)  are  analyzed,  with  re¬ 
sults  plotted  as  histograms  of  the  distribution  from  0-dB 
difference,  plus  plots  of  the  mean  of  that  distribution  and  the 
level  exceeded  by  90%  of  the  data.  Permits  user-selected 
constraints  similar  to  those  of  PLOT_DIS,  and  requires 
FORMIT  trip  files  and  predicted  signal  levels  as  inputs.. 

Several  data-reduction  program  changes  and  enhancements  were  implemented  for 
use  in  the  project.  These  changes  included: 

•  In  FORMIT,  input  was  modified  from  an  interactive  method  to  that 
of  a  constant  namelist  filename  (FORMIT.DAT). 

•  In  FORMIT,  automatic  extraction  of  ship’s  start  and  end  dates  and 
times  from  the  CNV  files  was  added. 
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•  In  PLOTJDIS  and  TRIPSTAT,  the  SUBDB  and  DELTDB  refer¬ 
ences  were  deleted. 

•  In  PLOT.DIS,  the  redundant  KFREQ  parameter  was  deleted  from 
namelist  and  the  function  was  automatically  generated  from  the 
FRQ  input.  Missing  ordinate  and  abscissa  labels  were  added. 

•  In  FORMIT,  PLOT.DIS  and  TRIPSTAT,  laser  printer  and  land¬ 
scape  switches  were  added. 

•  In  FORMIT,  PLOT.DIS  and  TRIPSTAT,  optional  output  spooling 
to  file  was  added,  permitting  later  plotting  on  a  laser  printer  in  the 
batch  mode,  with  print  enhancements  by  user  (llnewidth,  line  den¬ 
sity,  size,  x-y  position). 

•  In  PLOTJDIS,  capability  for  automatic  plotting  of  predicted  data 
from  a  FASTMC.ASC  file,  atop  recorded  data  plot,  was  added. 

•  In  TRIPSTAT,  capability  for  automatic  input  of  predicted  data  from 
a  FASTMC.ASC  file  was  added. 

•  As  initially  used,  PLOT.DIS  and  TRIPSTAT  permitted  ambiguous 
and  undetected  trip  number  errors.  For  example,  with  FORMIT 
data  files  for  trips  1-20  in  the  working  directory,  a  call  for  data  from 
trips  1-10  or  trips  11-20  yielded  plots  of  trip  1-10  data,  but  were 
marked  as  called.  Tht  programs  have  now  been  corrected. 

•  As  received,  PLOT.DIS  and  TRIPSTAT  were  capable  of  handling 
data  from  only  20  and  19  trips,  respectively.  This  limitation  has 
been  removed. 

•  PLOT.DIS  and  TRIPSTAT  have  been  compiled  for  32-bit  opera¬ 
tion,  speeding  execution  of  each  by  more  than  two-to-one  with 
respect  to  the  original  8-bit  compiles. 


13 


The  analysis  project  also  converted  all  TRIP.NRM  files  of  recorded  data  and  the  ac¬ 
companying  CNV.DAT  files  of  ship's  movement  from  the  52  numbered  trips,  and 
installed  the  files  on  project  computers.  Of  the  52  numbered  trips  of  the  GTS  Callaghan, 
data  from  trips  35  and  42  were  found  to  be  missing.  Therefore,  nighttime  analysis  was 
based  on  data  from  the  remaining  50  trips. 

A  listing  of  the  50  trips  of  the  GTS  Callaghan  is  given  in  table  1,  on  the  next  page. 
The  start  and  end  dates  and  times  (in  Universal  Time)  differ  slightly  from  those  listed 
in  table  1  of  the  Callaghan  report,  for  trips  1  through  22.  The  dates  and  times  of  table 
1  were  taken  directly  from  the  CNV  files  of  ship’s  movement,  and  were  utilized  in  the 
analyses  unless  recorded  data  near  the  origin  and/or  destination  were  found  to  be  in* 
appropriate. 
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Table  1.  GTS  Callaghan  Trips. 


TRIP 

ORIGIN 

DESTINATION 

START 

END 

i 

BAYONNB,  SC 

BREMMBRHAVEN 

03-39*66 

2230 

04-06-68 

1700 

2 

BREMERHAVEN,  W,  OR, 

CHARLESTON 

04-11.86 

2200 

04-19-86 

09Q0 

3 

''HARLESTON,  SC 

BAYONNB 

04-20-66 

0600 

04-21-86 

0700 

4 

BAYONNB 

BREMERHAVEN 

04-24-86 

0400 

08-02-66 

1018 

8 

BREMERHAVEN 

BAYONNE 

06-08-68 

0030 

06-13-86 

0830 

6 

BAYONNE 

BREMERHAVEN 

06-16-86 

0242 

08-32-66 

0600 

r 

BREMERHAVEN 

BAYONNE 

06-23-88 

2100 

06-30-68 

0600 

8 

BAYONNE 

BREMERHAVEN 

08-31-88 

2300 

06-07-88 

0830 

e 

BREMERHAVEN 

BAYONNB 

06-08-86 

2000 

06-14-86 

2000 

10 

BAYONNE 

BREMERHAVEN 

06-17-68 

0400 

06-33-88 

0800 

ii 

BHEMSRHAVBN 

BAYONNE 

06-28-86 

0400 

07-01-88 

0730 

13 

BAYONNE 

CHARLESTON 

07-03-86 

0100 

07-03-88 

0700 

13 

charleston 

SAVANNAH 

07-04-88 

0800 

07-04-86 

0930 

14 

BAVANNAH 

ALEXANDRIA,  EGYPT 

07-08-88 

0430 

07-19-68 

0230 

18 

ALEXANDRIA 

BREMERHAVEN 

07-21-66 

1400 

07-37-88 

2130 

16 

BREMERHAVEN 

BAYONNE 

07-29-66 

0000 

06-04-88 

1718 

IT 

BAYONNB 

BREMERHAVEN 

08-08-86 

0100 

08-14-6$ 

2200 

18 

BREMERHAVEN 

ALEXANDRIA 

08-16-88 

0200 

06-22-68 

0330 

19 

ALEXANDRIA 

SAVANNAH 

06-28-66 

0700 

00-06-68 

0712 

30 

SAVANNAH,  OA 

BREMERHAVEN 

09-10-66 

0330 

09-18-88 

1900 

21 

BREMERHAVEN 

NORFOLK 

09-20-66 

2400 

09-28-68 

2000 

23 

NORFOLK,  VA 

CHARLESTON 

09-30-88 

0119 

00-30-66 

1900 

23 

NORFOLK 

CHARLESTON 

10-33-86 

1600 

10-24-86 

0900 

24 

CHARLESTON 

ROTTERDAM,  HOLLAND 

10-36-68 

1300 

11-01-68 

1700 

28 

ROTTERDAM 

LIVERPOOL 

11-03-88 

1443 

11-04-86 

3000 

26 

LIVERPOOL,  ENGLAND 

ANTWERP 

11-06-88 

1834 

11-07-66 

1722 

27 

ANTWERP,  BELGIUM 

BREMERHAVEN 

11-08-86 

2400 

11-09-68 

0630 

28 

BREMERHAVEN 

CHARLESTON 

11-12-88 

0300 

11-21-88 

0830 

29 

CHARLESTON 

BAYONNB 

11-32-86 

0800 

11-23-88 

1100 

30 

BAYONNE 

NORFOLK 

11-28-88 

0130 

1 1  -28-68 

1000 

31 

NORFOLK 

SAVANNAH 

12-01-88 

2300 

12-02-68 

1830 

32 

SAVANNAH 

CRISTOBAL,  PANAMA 

12-06-88 

1430 

12-08-88 

0641 

33 

CRISTOBAL 

SAVANNAH 

12-10-88 

1740 

12-13-68 

0630 

34 

SAVANNAH 

PUERTO  CORTEZ,  HONDURAS 

12-17-88 

0200 

13-19-88 

1122 

36 

NO  DATA 

36 

CRISTOBAL 

JACKSONVILLE 

12-23-68 

2214 

12-20-68 

1100 

37 

JACKSONVILLE.  FL 

ZEBBRUQGE 

12-29-88 

0700 

01-06-88 

1646 

38 

2EEBAUQGB,  BELGIUM 

BREMERHAVEN 

01-08-66 

0vS4 

01-06-88 

1928 

39 

BREMERHAVEN 

BAYONNE 

01-10-66 

230C 

01-19-66 

1300 

40 

BAYONNE 

CHARLESTON 

02-04-66 

2330 

02-06-86 

1118 

41 

CHARLESTON 

ROTTERDAM 

02-08-66 

0400 

02-16-80 

0184 

43 

NO  DATA 

43 

BREMERHAVEN 

BEAUMONT 

02-22-66 

0600 

03-04-86 

0718 

44 

BEAUMONT,  TX 

CHARLESTON 

03-07-86 

1630 

03-06-66 

0918 

48 

CHARLESTON 

ROTTERDAM 

03-09-86 

0600 

03-16-86 

0800 

46 

ROTTERDAM 

BOGEN,  NORWAY 

03-16-66 

0300 

03-20-66 

1000 

47 

BOGEN 

MOREHEAD  CITY,  NC 

03-24-86 

0016 

04-02-86 

1300 

48 

MOREHEAD  CITY 

CHARLESTON 

04-08-86 

0016 

04-08-66 

0900 

49 

CHARLESTON 

BAYONNE 

04-06-86 

0800 

04-07-86 

0830 

80 

BAYONNE 

ROTTERDAM 

04-06-66 

0100 

04-18-86 

0900 

81 

ROTTERDAM 

BREMERHAVEN 

04-18-86 

0230 

04-16-66 

1101 

83 

BREMERHAVEN 

BAYONNE 

04-16-66 

0100 

04-28-66 

1700 
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SECTION  II-4 

LWPC  COMPUTER  PROGRAM 


LWPC,  Version  0,  as  written  for  the  VAX-VMS  system,  has  been  converted  for  use 
on  MS-DOS  computers  by  effort  in  this  project  and  in  associated  programs  for  the  De¬ 
fense  Nuclear  Agency,  LWPC  was  compiled  for  32-bit  operation  on  80386  machines 
using  the  Lahey  F-77L/EM-32  Fortran  compiler,  and  executes  at  speeds  comparable  to 
or  exceeding  those  of  the  DEC  MicroVAX*II  computer.  Specific  LWPC  program  mod¬ 
ules  available  for  project  use  are: 

•  PRESEG 

•  MODEFNDR 

•  SEGMWVGD 

•  MF.SW 

•  MERGE 

•  FASTMC 

•  PLTFMC 

LWPC  is  executed  on  MS-DOS  computers  using  the  VCL  (Virtual  Command  Lan¬ 
guage)  program,  which  creates  an  environment  similar  to  that  of  VAX  computers  using 
DCL  (DEC  Command  Language).  Thus,  the  various  DCL  command  files  created  by 
LWPC  in  its  automatic  sequences  may  be  used  and  interpreted  without  change.  Oper¬ 
ation  follows  descriptions  in  the  NOSC  LWPC  User's  Guide,  as  supplemented  by  a  con- 
tractor-written  User’s  Guide  for  MS-DOS  computers. 

Use  of  the  LWPC  programs  was  validated  in  two  ways.  First,  a  reference -path  cal¬ 
culation  was  produced  using  input  data  from  the  LWPC  User’s  Guide.  A  plot  from  this 
calculation  was  determined  to  be  equivalent  to  figure  D-3  of  the  NOSC  User's  Guide. 
Second,  the  LWPC  was  used  to  compute  signal  vs.  distance  data  for  the  four  curves  of 
VLF  predicted  data  shown  in  the  Callaghan  Report.  Computations  for  these  four  paths 
were  compared  with  those  of  the  Callaghan  report  and  documented  in  an  earlier  project 
report*. 


•"Comparative  Analysis  -  Callaghan  Data",  TCS  Report  TP-89-177,  31  May  1989. 
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A  number  of  LWPC  program  changes  and  enhancements  were  implemented  for  use 
in  the  project.  These  included: 

•  Automatic  operation  under  emulated  VAX-VMS  environment. 

•  Filename  conventions  revised  to  accommodate  the  8  +  3  character 
filename  limitation  of  MS-DOS. 

•  Automatic  generation  of  Norton  FILEINFO.FI  files  to  annotate  re¬ 
stricted  filenames.  Initial  implementation  adds  path  bearings  and 
date/time  tags. 

•  In  PRESF.G,  four-digit  power  levels  may  now  be  specified.. 

•  In  PRESEG,  the  NPRINT  function  war-  added  for  automatic  feed 
to  the  FASTMC  command  file, 

•  In  FASTMC  and  PLTFMC,  Amplitude/Phase  vs.  Distance  print  for¬ 
mat  (for  NPRINT  >0)  was  changed  to  nine  columns  (instead  of 
twelve)  for  easier  printing  on  80-column  printers. 

•  In  FASTMC,  a  two-column  listing  of  distance  and  amplitude  in  an 
automatically-generated  ,ASC  file  was  added,  for  use  in  inputting 
predicted  data  automatically  to  PLOT_DIS  and  to  TRIPSTAT. 

•  In  PLTFMC,  NRPTS  was  added  to  namelist  input  parameters. 

•  In  PLTFMC,  landscape  mode  was  added. 

•  In  PLTFMC,  an  output  switch  was  added  to  permit  selection  of  plot¬ 
ter  or  laser  printer  HPGL  emulation. 

•  In  PLTFMC,  output  spooling  to  file  was  added,  permitting  later 
plotting  on  a  laser  printer  in  batch  mode,  with  print  enhancements 
by  the  user  (linewidth,  line  density,  size,  x-y  position,  etc.). 


•  Automatic  configuration  management  and  archiving  of  all  program 
versions  by  number  was  implemented. 
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SECTION  11*5 

MEASURED  NIGHTTIME  DATA 
II-5.1  SELECTED  TRIP  SETS 

Because  of  the  relative  variability  of  recorded  nighttime  signal  data,  with  respect  to 
the  more  uniform  signals  of  the  daytime  environment,  it  was  important  to  constrain  the 
Callaghan  measurements  as  much  as  possible  before  comparing  with  the  LWPC  predic¬ 
tions.  Normalized  data  from  each  trip  was  examined  for  insured  integrity  at  each 
frequency  of  interest,  rejecting  those  cases  where  obvious  anomolies  in  the  measure¬ 
ments  had  occurred  or  where  the  ship’s  course  was  outside  the  area  represented  by  the 
selected  propagation  path  to  be  computed  by  the  prediction  model.  The  data  reduction 
program  PLOT_DIS  was  run  for  individual  trip  data  at  each  frequency,  and  a  reference 
manual  of  plotted  data  was  assembled  as  an  aid  to  selecting  the  groups  of  valid  trips. 

Individual  plots  of  the  ship’s  course  for  each  trip  were  also  made,  using  the  Callaghan 
navigation  files.  The  NOSC  program  WMAP  was  modified  to  permit  automatic  plotting 
of  the  ship’s  course  on  a  world  map,  and  a  reference  manual  of  individual  trip  plots  was 
created.  (Figure  2  is  an  example  of  such  plots  for  a  group  of  trips.) 

A  “Primary  Trip  Set”  of  recorded  data  was  then  selected,  following  the  procedure 
diagramed  in  figure  4,  and  repeating  the  process  for  each  of  the  four  frequencies.  Trip 
data  files  containing  other  than  all-uighttime  data  were  first  eliminated.  Plots  of  the 
ship’s  course  for  each  trip  were  then  examined,  and  trips  out  of  the  North  Atlantic  area 
of  principal  interest  were  rejected,  along  with  several  short  trips  in  the  European  area 
and  between  CONUS  ports.  Some  of  the  trips  to  Central  America  were  retained  for 
the  16.0-  and  19.0-kHz  British  transmitters,  since  such  trips  appeared  to  provide  useful 
data  in  the  6-  to  8-Mm  range.  The  Central  America  trips  were  rejected  for  the  CONUS- 
based  transmitters,  however,  since  such  locations  were  considerably  removed  from  the 
propagation  paths  of  interest.  Finally,  the  remaining  trip  data  was  examined  for  integ¬ 
rity,  rejecting  files  that  were  obviously  contaminated  by  noise. 

The  remaining  data  files  were  defined  as  the  Primary  Trip  Sets,  and  aggregated  re¬ 
corded  data  from  these  sets  were  used  for  the  comparisons.  Tables  2-5  on  the  following 
pages  record  the  final  trip  selections,  and  indicate  the  reasons  for  rejecting  individual 
data  files  using  the  reject  codes  of  figure  4. 
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Figure  4.  Methodolog>  diagram,  Callaghan  trip  re  corded  data, 
selection  of  measured  data  for  comparison. 
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Table  2.  Primary  trip  set,  16.0  kHz  (large  type) 


REJECT 

TRIP 

ORIGIN 

DESTINATION 

START 

END 

CODE 

i 

BAYONNE,  SC 

BREMMERHAVEN 

03-29-89 

2230 

0**0M5 

1700 

2 

BREMERHAVEN,  W,  OR, 

CHARLESTON 

04-11-68 

2200 

04-19-66 

0900 

2 

3 

CHARLESTON,  1C 

BAYONNE 

042010 

0400 

04-31 -IS 

0700 

A 

BAYONNE 

BREMERHAVEN 

04-24-99 

0400 

05-02-89 

1018 

9 

BREMERHAVEN 

BAYONNE 

09-05-89 

0030 

08-13-69 

0830 

4 

fl 

IAYONNI 

BREMERHAVEN 

08-  15-M 

0243 

ami 

0400 

4 

7 

SRSMERHAVEN 

BAYONNE 

0B-23-AB 

2100 

M-30-1S 

0100 

4 

i 

BAYONNE 

BREMERHAVEN 

SS-H-IB 

MOO 

0S-O7-AS 

0130 

e 

BREMERHAVEN 

BAYONNE 

09-08-89 

2000 

06-14-69 

2000 

10 

BAYONNE 

BREMERHAVEN 

08-17.68 

0400 

06-23-68 

0600 

ii 

bREMERHAVEN 

BAYONNE 

09-29-68 

0400 

07-01-68 

0730 

2 

•t 

EA  TONNE 

CHARLIETON 

otoe-ii 

0130 

07.03-S! 

0700 

1 

M 

CHARLEITON 

SAVANNAH 

97-04-11 

0603 

07-04 -SB 

0430 

14 

SAVANNAH 

ALEXANDRIA  EGYPT 

07-06-66 

0430 

07-19-66 

0230 

19 

ALEXANDRIA 

BREMERHAVEN 

07-21-69 

1400 

07-27-88 

2130 

16 

BREMERHAVEN 

BAYONNE 

07-29-69 

OOOO 

08-04-6B 

1718 

17 

BAYONNE 

BREMERHAVEN 

06-06-69 

0100 

06-14-69 

2200 

2 

'1 

■AIMERHAVEN 

ALEXANDRIA 

01-11-19 

0100 

01-U-SS 

0330 

19 

ALEXANDRIA 

8AVANNAH 

08-29-69 

0700 

09-06-89 

0712 

20 

SAVANNAH.  OA 

BREMERHAVEN 

09-10-66 

0330 

09-16-69 

1900 

21 

BREMERHAVEN 

NORFOLK 

09-20-68 

2400 

09-26-66 

2000 

2 

» 

NOR'OLK.  VA 

CHARLESTON 

0*3018 

0114 

04-30-19 

'100 

2 

1) 

NORFOLK 

CHARLESTON 

10-13-19 

1400 

10-14-11 

0100 

24 

CHARLESTON 

ROTTERDAM,  HOLLAND 

10-26-86 

1300 

11-01-86 

1700 

2 

it 

ROTTERDAM 

LIVERPOOL 

11-03-11 

1444 

11-04-45 

MOO 

2 

» 

UVBRROOL,  1 NO LAND 

ANTWERP 

11-OS-SI 

1S34 

11-07-16 

1TM 

2 

r 

ANTWIRR,  aiLOlUM 

BREMERHAVEN 

11-OS-IS 

1400 

11*04*11 

ouo 

26 

BREMERHAVEN 

CHARLESTON 

11-12-66 

0300 

11-21-86 

0930 

2 

n 

CHARLEITON 

BAYONNE 

ii-n-sB 

osoo 

1 1 -13-11 

1100 

2 

JO 

BAYONNE 

NORFOLK 

11-IS-SS 

0100 

'  1-14-41 

1000 

1 

31 

NORROLX 

IAVANNAM 

1S-01-SS 

3300 

1301-11 

1110 

32 

SAVANNAH 

CRISTOBAL  PANAMA 

12-06-86 

1430 

12-06*85 

0641 

33 

CRISTOBAL 

SAVANNAH 

12-10-08 

1740 

12-13-66 

0630 

34 

SAVANNAH 

PUERTO  CORTEZ  HONDURAS 

12-17-B5 

0200 

12-19-85 

1122 

1 

33 

38 

NO  OATA 

CRISTOBAL 

JACKSONVILLE 

12-23*05 

2214 

12-29-85 

1100 

37 

JACKSONVILLE,  PL 

2EEBRU0QE 

12-29-85 

0700 

01-08-86 

1846 

1 

31 

ZSIBRUOOE,  IILQIUM 

BREMERHAVEN 

01-04*44 

0454 

01-04*41 

ms 

39 

BREMERHAVEN 

BAYONNE 

01-1048 

2300 

01-19-66 

1300 

2 

*0 

BAYONNE 

CHARLEITON 

03-04-44 

4330 

01-04-M 

111B 

1 

41 

charuiiton 

ROTTERDAM 

03-04-44 

0400 

01-1MB 

01M 

1 

a 

43 

NO  DATA 

BREMERHAVEN 

BEAUMONT 

02-22-66 

0600 

03-04-66 

0719 

2 

44 

BEAUMONT,  TX 

CHARLEITON 

03-0746 

1100 

03-OM1 

0411 

48 

CHARI  ESTON 

ROTTERDAM 

03-09-66 

0600 

03-16-66 

0800 

2 

*1 

ROTTERDAM 

EOaEN,  NORWAY 

03-11-44 

0300 

30-n-ss 

1000 

47 

BOOEN 

MOREHEAD  CITY,  NC 

03-24-66 

0018 

04-02-86 

1300 

2 

<1 

MOREhBAD  CITY 

CHARLEITON 

0448-46 

001S 

34-01-M 

0100 

2 

*• 

CHARLESTON 

IAYONNI 

04*0v*«i 

0900 

34-07-41 

0110 

60 

BAYONNE 

ROTTERDAM 

04-06*88 

0100 

04-15-86 

0900 

1 

3' 

ROTTERDAM 

WIMINHAVIN 

04-14-44 

0130 

04-11-11 

1101 

92 

BREMERHA'  EN 

BAYONNE 

04-16-66 

0100 

04-2640 

1700 

20 


Table  3.  Primary  trip  set,  19.0  kHz  (large  type). 


REJECT 

CODE 

TRIP 

ORIGIN 

DESTINATION 

START 

END 

i 

8AYONN5,  SC 

BREMMERHAVEN 

03.28.88 

,2230 

04-05-88 

1700 

2 

BREMERHAVEN,  W.  QR, 

CHARLESTON 

04.11-85 

2200 

04-18-89 

0800 

2 

a 

0H4BUIT0N.  1C 

•AYONNI 

ok-m-m 

0900 

04-11-91 

0700 

4 

4 

■AYONNI 

■MUIAhAVIN 

04-84-90 

0400 

oi-oa-ii 

1014 

9 

BREMERHAVEN 

BAYONNE 

06-08-88 

0030 

05-13-85 

0830 

6 

BAYONNE 

BREMERHAVEN 

06- 15-85 

0242 

09-22-85 

0800 

r 

BREMERHAVEN 

BAYONNE 

05-23-85 

2100 

05-30-85 

0800 

8 

BAYONNE 

BREMERHAVEN 

05-31-65 

2300 

08-07-85 

0630 

0 

BREMERHAVEN 

BAYONNE 

09-08-86 

2000 

00-14-89 

2000 

10 

BAYONNE 

BREMERHAVEN 

08-17-85 

0400 

08-23-86 

0800 

11 

BREMERHAVEN 

BAYONNE 

05-25-85 

0400 

07-01-68 

0730 

2 

It 

•AYONNI 

CHAAlUTON 

QT-et-ii 

0100 

07-09-19 

OTOO 

1 

*3 

CHANJITON 

•AVANNAH 

07-04*19 

0900 

07-04-99 

0990 

14 

SAVANNAH 

ALEXANDRIA,  EGYPT 

07-08-88 

0430 

07-18-65 

0230 

19 

ALEXANDRIA 

BREMERHAVEN 

07-21-88 

1400 

07-27-88 

2130 

19 

BREMERHAVEN 

BAYONNE 

07-28-86 

0000 

08-04-88 

1718 

17 

BAYONNE 

BREMERHAVEN 

06-06-66 

0100 

06-14-16 

2200 

2 

14 

■MMIRHAVIN 

AltXANOAIA' 

09-1  Ml 

OHO 

01-lt-tt 

0990 

10 

ALEXANDRIA 

SAVANNAH 

08-28-68 

0700 

08-06-18 

0712 

20 

SAVANNAH,  OA 

BREMERHAVEN 

08-10-66 

0330 

08-18-98 

1800 

21 

BREMERHAVEN 

NORFOLK 

08-20-88 

2400 

08-26-68 

2000 

2 

It 

NOMOLK.  VA 

CHAALMTON 

OI-IO-II 

0119 

09-90-19 

1900 

2 

13 

NOMOLK 

CHAAUITON 

10-13-11 

1900 

10-14-91 

0900 

24 

CHARLESTON 

ROTTERDAM,  HOLLAND 

10-26-88 

1300 

11-0148 

1700 

2 

II 

AOTTinOAM 

UVIAPOOI. 

11-09-11 

1441 

11-04-99 

1000 

2 

If 

UVM40CL  1  NOLAND 

ANTWtM 

1 1  -04-91 

1194 

11-07-91 

1711 

2 

ST 

ANTWIW,  IILUIUM 

•AIMIAHAVIN 

11-OMJ 

1400 

11-09*99 

0990 

28 

BREMERHAVEN 

CHARLESTON 

11-12-86 

0300 

11-2148 

0B30 

2 

» 

CHAAlWON 

•AYONNI 

n-lt-ti 

OIOO 

11-19-99 

1100 

2 

ae 

■AYONNI 

NOMOLK 

1 1 -till 

0190 

11-19-99 

1000 

1 

91 

NOMOLK 

IAVANNAH 

11-01-11 

1900 

H-0I-9B 

1490 

32 

SAVANNAH 

CRISTOBAL.  PANAMA 

12-05-85 

1430 

12-0648 

0641 

33 

CRISTOBAL 

SAVANNAH 

12-10-85 

1740 

12-13-85 

0630 

34 

SAVANNAH 

PUERTO  CORTEZ,  HONDURAS 

12-17-85 

0200 

12-18-85 

1122 

1 

at 

NO  DATA 

39 

CRISTOBAL 

JACKSONVILLE 

12-23-65 

2214 

12-264S 

1100 

37 

JACKSONVILLE,  PL 

ZEEBRUOQE 

12-28-85 

0700 

01-08-66 

1648 

1 

31 

ziimuaai,  iiloium 

■A1M1RHAVIN 

01-09-49 

0904 

01-09-49 

1M4 

30 

BREMERHAVEN 

BAYONNE 

01-10-86 

2300 

01-18-80 

1300 

2 

40 

•AYONNI 

OMAALIITON 

0I-O4-M 

1930 

0903-94 

1114 

1 

41 

CHAAlJITON 

AOTTISOAM 

0109-91 

0400 

09-14-99 

01*4 

1 

41 

NO  OATA 

43 

BREMERHAVEN 

BEAUMONT 

02-22-66 

0800 

03-04-86 

0718 

2 

44 

IIAUMONT,  TO 

ohaaliiton 

03-0T-II 

1190 

0904-14 

0«4 

48 

CHARLESTON 

ROTTERDAM 

03-08-88 

0800 

03-18-88 

0600 

2 

41 

AOTT«AOAM 

I09IN,  NOWAY 

03-11-91 

0300 

09-10-14 

1000 

47 

BOOEN 

MOREHEAD  CITY,  NC 

03-24-88 

0015 

04-02-86 

1300 

2 

41 

MOAIHIAD  CITY 

CHANLllTON 

04-01-99 

0019 

94-09-14 

0K0 

2 

41 

CHAAlMTON 

•AYONNI 

04-09-99 

0900 

94-07-14 

0410 

90 

BAYONNE 

ROTTERDAM 

04-08-59 

0100 

04-1646 

0800 

1 

n 

AOTTtAOAM 

■MMINhAVIN 

04-19-99 

0190 

04-14-99 

1101 

92 

BREMERHAVEN 

BAYONNE 

04-18-86 

0100 

04-28-88 

1700 

21 


Table  4.  Primary  trip  set,  21.4  kHz  (large  type). 


REJECT 

TRIP 

ORIGIN 

DESTINATION 

START 

END 

CODE 

i 

BAYONNE,  SC 

BREMMERHAVEN 

03-29-68 

2230 

04*06*88 

1700 

2 

BREMERHAVEN,  W,  OR, 

CHARLESTON 

04*  11.  as 

2200 

04*18*08 

0800 

2 

3 

CHA4LIIT0N,  ic 

IAVONNI 

0410  13 

0400 

04-11-31 

0700 

4 

BAYONNE 

BREMERHAVEN 

04-24.88 

0400 

08-02*88 

1018 

S 

BREMERHAVEN 

BAYONNE 

06-08-88 

0030 

05-13-88 

0630 

6 

BAYONNE 

BREMERHAVEN 

08.18-66 

0242 

08-22-88 

0800 

7 

BREMERHAVEN 

BAYONNE 

08-23-88 

2100 

08-30-66 

0600 

a 

BAYONNE 

BREMERHAVEN 

08.31*88 

2300 

06*07-86 

0630 

6 

BREMERHAVEN 

BAYONNE 

08-08-88 

2000 

06*14*88 

2000 

10 

BAYONNE 

BREMERHAVEN 

06*  17*88 

0400 

06*23*66 

0800 

11 

BREMERHAVEN 

BAYONNE 

06*28  88 

0400 

07*01*88 

0730 

2 

•a 

IAVONNI 

CHARLIITON 

07  01-33 

0100 

07-03-30 

0700 

1 

'3 

ChABLIITON 

IAVANNAH 

07*04.18 

0400 

07-04-31 

0130 

14 

SAVANNAH 

ALEXANDRIA  EQYRT 

07*0846 

0430 

07-16-18 

0210 

ie 

ALEXANDRIA 

BREMERHAVEN 

07.21*88 

1400 

07*27*86 

2130 

10 

BREMERHAVEN 

BAYONNE 

07.28*88 

0000 

08*04*86 

1716 

17 

BAYONNE 

BREMERHAVEN 

08-06-88 

0100 

06-14*86 

2200 

2 

ii 

■MMWHAVIN 

ALIXAN0MA 

M-1I-U 

ono 

0A-lt.ll 

0Q30 

18 

ALEXANDRIA 

SAVANNAH 

08-28-88 

0700 

08*06-86 

0712 

20 

SAVANNAH,  OA 

BREMERHAVEN 

08*1046 

0330 

09*1848 

1900 

21 

BREMERHAVEN 

NORFOLK 

08*2046 

2400 

09*21*16 

2000 

2 

at 

NOMOLK,  VA 

OMAAlHTON 

Ot-lO-ll 

out 

09*Q0*II 

two 

2 

13 

NOBBOLK 

CHARLIITON 

101311 

1400 

10-14-11 

0400 

24 

CHARLESTON 

ROTTERDAM,  HOLLAND 

10*2846 

1300 

11-0146 

1700 

2 

II 

AommoAM 

UVIWOOl 

I1-«*M 

1441 

11-04-11 

MOO 

2 

H 

UVIMBOOL,  INQIANO 

ANTWIM 

M -04-31 

1U4 

1 1*47*31 

ITU 

2 

IT 

ANTWIHB,  IILQIUM 

IHIMINAVIN 

11-OMI 

1400 

11-04-tl 

0430 

21 

BREMERHAVEN 

CHARLESTON 

11-1248 

0300 

11*2146 

0830 

2 

H 

OHAAUITON 

IAVONNI 

1 1  'li*il 

0400 

11 -11-31 

1100 

2 

10 

IAVONNI 

NOMOLK 

1  MI-11 

01 M 

1 1  *11-31 

1000 

1 

1* 

NOMOLK 

IAVANNAH 

11-OMI 

ttoo 

it-OI-ll 

1130 

3 

11 

iavannah 

CBiITOIAL  RANAMA 

11-03  IB 

1410 

11-01-10 

0141 

3 

33 

CBiITOIAL 

IAVANNAH 

11-10-41 

<T40 

11-11-30 

0410 

3 

14 

IAVANNAH 

AUINTO  COKTU,  “CNOURAI 

1I-1T-4I 

otoo 

11-11-10 

tin 

1 

11 

NO  DATA 

3 

M 

CBIITOIAL 

JACKIONVILU 

11-13-41 

1114 

U-I4-II 

1100 

37 

JACKSONVILLE,  RL 

ZEEBRUQQE 

12-2848 

0700 

01*06*86 

1846 

1 

M 

ZIIIBU90I,  IILQIUM 

IBIMIBHAVIN 

01-04-44 

0444 

01-00-34 

1MI 

38 

BREMERHAVEN 

BAYONNE 

01*10*86 

2300 

01*18*86 

1300 

2 

40 

■AVONNI 

CHARLIITON 

0104-44 

two 

94-04  34 

1114 

1 

41 

OHAHH8TON 

noniAOAM 

010404 

04C0 

01-1434 

0114 

1 

41 

NO  DATA 

43 

BREMERHAVEN 

BEAUMONT 

02-22-88 

0600 

03*04-86 

0718 

2 

44 

IIAUMONT,  TX 

CHARLIITON 

0107*44 

1130 

03-64-34 

0414 

48 

CHARLESTON 

ROTTERDAM 

03-0948 

0800 

03*16*86 

0600 

2 

4« 

BOTTIBDAM 

IOQIN,  NORWAV 

01-11-44 

OKO 

03-10-14 

1000 

47 

BOQEN 

MOREHEAD  CITY,  NC 

03*24-06 

0018 

04*02*66 

1300 

2 

41 

momhiao  orrv 

CHARLIITON 

0404-14 

OO'I 

04.00-34 

OMO 

2 

II 

ChARLIITON 

•AVONNI 

D4-0I19 

0100 

04-07-14 

0410 

80 

BAYONNE 

ROTTERDAM 

04*06-66 

O’CO 

04-1646 

0800 

1 

it 

ROTTfROAM 

IMMI4HAVIN 

04-1404 

0133 

04-14-14 

1101 

82 

BREMERHAVEN 

BAYONNE 

04*16*66 

0100 

04-25*86 

1700 

22 


Table  5.  Primary  trip  set,  24.0  kHz  (large  type). 


REJECT 

CODE 

TRIP 

ORIGIN 

DESTINATION 

START 

END 

1 

BAYONNE,  SC 

BREMMERHAVEN 

03-29-98 

2330 

04-0845 

1700 

a 

BREMENHAVEN,  W.  OR, 

CHARLESTON 

04-11-85 

2200 

04-1048 

0800 

a 

a 

OHARUITON,  90 

IAYONNI 

044049 

0900 

044149 

0700 

4 

BAYONNE 

BREMERHAVEN 

04-44-85 

0400 

06-02-85 

1018 

S 

BREMSRHAVEN 

BAYONNE 

05-05-85 

0030 

05-1349 

0830 

9 

BAYONNE 

BREMERHAVEN 

08-15-85 

0242 

05-22-85 

0800 

7 

BREMERHAVEN 

BAYONNE 

05-23-85 

2100 

05-30-86 

0800 

6 

BAYONNE 

BREMERHAVEN 

05-31-85 

2300 

09-0745 

0630 

9 

BREMENHAVEN 

BAYONNE 

06-38-85 

2000 

08-1446 

2000 

10 

BAYONNE 

BREMERHAVEN 

OB-17-85 

0400 

08-2345 

0800 

11 

BREMENHAVEN 

BAYONNE 

05-25-85 

0400 

07-0149 

0730 

a 

ii 

IAYONNI 

OHARLIITON 

97-91-91 

0100 

0749-99 

9700 

i 

ia 

OHARUITON 

SAVANNAH 

97-04-91 

0900 

07-04-11 

0990 

14 

SAVANNAH 

ALEXANDRIA,  EQYPT 

07-08-85 

0430 

07-1845 

0230 

16 

ALEXANDRIA 

BREMERHAVEN 

07-2145 

1400 

07-2745 

2130 

19 

BREMENHAVEN 

BAYONNE 

07-29-85 

0000 

08*0445 

1715 

17 

BAYONNE 

BREMERHAVEN 

08-08-86 

0100 

08-14-85 

2200 

a 

il 

IRIMIRHAVIN 

ALIXANORIA 

09-19-99 

0M0 

09-9149 

0990 

19 

ALEXANDRIA 

SAVANNAH 

08-25-85 

0700 

00-0646 

0712 

ao 

SAVANNAH,  OA 

BREMERHAVEN 

09-1045 

0330 

08-1846 

1800 

ai 

BREMERHAVEN 

NORFOLK 

09-20-85 

2400 

08-2849 

2000 

a 

n 

NORFOLK,  VA 

OHARLIITON 

09*90*99 

0119 

09-90-11 

1900 

a 

n 

NORFOLK 

OHARLIITON 

10-19-99 

WO 

10-14-99 

0900 

24 

charleston 

ROTTERDAM,  HOLLAND 

10-8545 

1300 

11-0145 

1700 

a 

» 

ROTTIROAM 

UVIRROOL 

11 -03-91 

1441 

1144-99 

MOO 

a 

H 

UVIRROOU  INOLAND 

ANTWIRR 

11-09-99 

1934 

114749 

ITU 

a 

*7 

ANTWIRR,  IILOIUM 

■RIMIRHAVIN 

11-09-99 

1400 

1149-91 

0930 

as 

BREMERHAVEN 

CHARLESTON 

11-1848 

0300 

11-2145 

0830 

a 

M 

oharliiton 

IAYONNI 

11-1949 

0900 

11-1349 

1100 

a 

M 

BAYONNI 

NORFOLK 

11 -19-91 

0190 

11-1949 

1000 

i 

91 

NORFOLK 

SAVANNAH 

11-01-99 

1900 

114149 

1900 

3 

91 

SAVANNAH 

CRISTOBAL.  RANAMA 

11-09-91 

1490 

114941 

0941 

3 

99 

CRIITOIAL 

SAVANNAH 

11-1049 

1740 

19-1949 

0930 

3 

90 

SAVANNAH 

RUIRTO  OORTTB,  HONDURAS 

19-17-99 

0100 

11-1949 

1191 

1 

99 

NO  DATA 

3 

99 

CRISTOBAL 

JACKSONVUI 

19-1949 

ni4 

19-1949 

1100 

37 

JACKSONVILLE,  FL 

ZEEBRUQQE 

12-2845 

0700 

01-0948 

1648 

1 

99 

ZIISRU93I  IIL3IUM 

BRIMIRHAVIN 

01-09-99 

0994 

314949 

1999 

39 

BREMERHAVEN 

BAYONNE 

01-1049 

2300 

01-18-88 

1300 

2 

40 

BAYONNI 

OHARLIITON 

01-0449 

1930 

014949 

1119 

1 

41 

OHARLBSTON 

ROTTIROAM 

01-0949 

0400 

09-1949 

0194 

1 

49 

NO  DATA 

43 

BREMERHAVEN 

BEAUMONT 

03-2349 

0800 

03-0446 

0718 

a 

44 

BIAUMQNT,  TX 

OHARLIITON 

03-07-99 

1990 

034949 

0119 

48 

CHARLESTON 

ROTTERDAM 

03-0048 

0800 

03-16-88 

0800 

2 

49 

ROTTIROAM 

■OOIN,  NORWAY 

03-1949 

0900 

03-1049 

1300 

47 

BOQEN 

MOREHEAO  CITY,  NC 

03-2448 

0018 

04-02-88 

1300 

2 

49 

MORSHIAO  CITY 

OHARLIITON 

04-01-99 

nil 

0449-99 

0900 

2 

49 

OHARLIITON 

BAYONNI 

04,0949 

0900 

0447-99 

0930 

90 

BAYONNE 

ROTTERDAM 

04-08-86 

0100 

04-1545 

0800 

1 

91 

ROTTIROAM 

BRIMIRHAVIN 

04-19-99 

0130 

04-19-99 

1101 

92 

BREMERHAVEN 

BAYONNE 

04-1849 

0100 

04-25-88 

1700 

23 
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11-5,2  SEASONAL  EFFECTS 

An  attempt  was  made  to  separate  the  recorded  data  into  groups,  sorted  by  season. 
Table  6  shows  the  50  Callaghan  trips  divided  into  the  three-month  seasons  of  spring, 
summer,  autumn  and  winter,  with  the  four  periods  defined  by  the  summer  solstice  on 
June  21  and  the  winter  solstice  on  December  22.  Similarly,  six-month  “warm"  and 
“cold"  seasons  were  defined  around  the  same  dates. 

It  was  decided,  however,  that  the  Callaghan  database  was  too  small  to  support  mean¬ 
ingful  separation  by  season.  Figures  5-8  illustrate  the  recorded  data  from  the  Rugby 
transmitter  broadcasting  at  16  kHz,  separated  into  spring,  summer,  autumn  and  winter 
groupings,  respectively,  as  an  example  case.  Overlays  of  these  plots  revealed  similar  pat¬ 
terns  of  recorded  signal  with  distance -there  were  no  pronounced  seasonal  effects  that 
would  justify  the  use  of  different  ionospheric  models  in  the  LWPC. 

Similar  conclusions  were  reached  in  comparing  the  six-month  seasonal  data.  Again, 
the  great  variability  of  nighttime  data  due  to  undefined  changes  in  the  propagation  chan¬ 
nel  seemed  to  obscure  any  effects  that  were  purely  seasonal. 
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WARM  !  COLD  SEASON  i  WARM  SEASON 


Table  6.  Callaghan  trips  divided  by  season. 


TRIP 

ORIGIN 

DESTINATION 

START 

END 

i 

BAYONNE,  SC 

BRSMMERHAVEN 

03*29-99 

2230 

04*00*68 

1700 

2 

BREMERHAVEN,  W.  OR, 

CHARLESTON 

04.11-86 

2200 

04*19*68 

0900 

3 

charleston,  sc 

BAYONNE 

04.20*68 

0900 

04*21*88 

0700 

4 

BAYONNE 

BREMERHAVEN 

04-24.89 

0400 

08*02-68 

1016 

5 

BPEMERhAVEN 

BAYONNE 

08*06*88 

0030 

06*13*68 

0630 

e 

BAYONNE 

BREMERHAVEN 

08*18*68 

0242 

09-22-88 

0600 

7 

BPEMERHAVEN 

BAYONNE 

09-23-69 

2100 

08-30*89 

0600 

8 

BAYONNE 

BREMERHAVEN 

08*31*08 

2300 

06*07.08 

0830 

e 

BREMERHAVEN 

BAYONNE 

00*06*88 

2000 

00*14*66 

2000 

10 

BAYONNE 

BREMERHAVEN 

06*17*68 

0400 

06*23*68 

0600 

*  1 

BPEMERHAVEN 

BAYONNE 

06*28*68 

0400 

07-01  *68 

0730 

i: 

BAYONNE 

charleston 

07*02*68 

0100 

07.03*08 

0700 

13 

CHARLESTON 

savannah 

07*04*68 

0800 

07*04*88 

0930 

u 

SAVANNAH 

ALEXANDRIA,  B3YPT 

07*06*88 

0430 

07*19*68 

0230 

13 

ALEXANDRIA 

BREMERHAVEN 

07*21*69 

1400 

07*27*08 

2130 

•3 

BREMERHAVEN 

BAYONNE 

07*29*88 

0000 

09*04*98 

1710 

17 

baySnne 

brTmerhaven 

06*08*69 

0100 

06*14*88 

2200 

16 

BPEMERHAVEN 

ALEXANDRIA 

06*10*68 

0200 

06*32*86 

0330 

10 

ALEXANDRIA 

SAVANNAH 

06*28*68 

0700 

09*00*06 

0712 

20 

SAVANNAH,  OA 

BREMERHAVEN 

06*10*88 

0330 

09*18*86 

1900 

21 

Vr'emYrhaven" 

n6W5lk 

09*20-68 

2400 

09*26-88 

2000 

22 

NORFOLK,  VA 

CHARLESTON 

09*30*08 

0119 

09*30*68 

1900 

23 

NORFOLK 

CHARLESTON 

10*23*68 

1600 

10*24*98 

0600 

3* 

B»r*n 

mi'Un 

■uvl 

26 

ROTTERDAM 

LIVERPOOL 

11*03*68 

1448 

11*04*88 

2000 

20 

LIVERPOOL,  ENGLAND 

ANTWERP 

11*00*69 

1834 

11*07*08 

1722 

27 

ANTWERP,  BELGIUM 

BREMERHAVEN 

11*08*68 

2400 

11-09*08 

0630 

28 

BREMERHAVEN 

CHARLESTON 

11*12*88 

0300 

11*21*88 

0630 

20 

charleston 

BAYONNE 

11*22*66 

0600 

11*23*88 

1100 

30 

BAYONNE 

NORFOLK 

11*26*08 

0130 

11*20*68 

1000 

31 

NORFOLK 

SAVANNAH 

12*01*68 

2300 

12*02*88 

1630 

32 

SAVANNAH 

CRISTOBAL,  PANAMA 

12*06*88 

1430 

12*00*66 

0041 

33 

CRISTOBAL 

SAVANNAH 

12*10*09 

1740 

12*13*96 

0030 

34 

SAVANNAH 

PUERTO  CORTEZ,  HONDURAS 

12*17*68 

0200 

12*16*88 

1122 

39 

NO  DATA 

30 

CRISTOBAL 

JACKSONVILLE 

12*23*68 

2214 

12*20*88 

1100 

37 

JACKSONVILLE,  FL 

ZEEBRUQGE 

12*20*66 

0700 

01*00*80 

1646 

38 

2SEBRUQQE,  BELGIUM 

BREMERHAVEN 

01*06*60 

0084 

01*00*80 

1926 

30 

BREMERHAVEN 

BAYONNE 

01*10*60 

2300 

01*19*90 

1300 

40 

—  i  ti  if  i  ii  — — 

■H-M-gT-M 

1118 

41 

CHARLESTON 

ROTTERDAM 

02*06*60 

0400 

02*16*86 

0184 

42 

NO  DATA 

43 

BREMERHAVEN 

BEAUMONT 

03*22*60 

0000 

03*04*80 

0710 

44 

BEAUMONT,  TX 

CHARLESTON 

03*07*00 

1830 

03*00*80 

0916 

49 

charleston 

ROTTERDAM 

03*09*80 

0600 

03*16*60 

0800 

40 

’ROTTERDAM 

BOQEN,  NORWAY 

63*16*06 

0300 

03-20*86 

1000 

47 

BOQEN 

MOREHEAD  CITY,  NC 

03*24*08 

0016 

04*02*80 

1300 

49 

MORIHEAD  CITY 

CHARLESTON 

04*06*60 

0016 

04.06*80 

0900 

40 

charleston 

BAYONNE 

04*00*06 

0600 

04*07*80 

0630 

90 

BAYONNE 

ROTTERDAM 

04*00*60 

0100 

04*18*60 

0900 

91 

ROTTERDAM 

BREMERHAVEN 

04*10*80 

0230 

04*16*86 

1101 

92 

BREMERHAVEN 

BAYONNE 

04*18*69 

0100 

04*28*86 

1700 

25 


SPRING  WINTER  AUTUMN  SUMMER  SPRING 


TCS  PLOT _PIS  09/26/89  16:  14 

Freq:  16.000  kfreq:  160 

NIGHT  (29- Mar-85)  to  (25-Apr-B6)  deldbl:  B.  deldb2:  12. 


Figure  5.  Typical  recorded  data  for  Spring  season,  Rugby  at  16.0  kHz 
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Figure  6.  Typical  recorded  data  for  Summer  season,  Rugby  at  16.0  kHz. 


II-5.3  SOLAR  EVENT  ANALYSIS 

It  was  recognized  that  unusual  solar  activity  could  have  occurred  during  the  Callag¬ 
han  trips,  resulting  in  ionospheric  and/or  geomagnetic  field  disturbances  and  consequent 
anomolous  signal  strength  recordings,  To  investigate  this  possibility,  solar  activity  re¬ 
cords  were  obtained  from  the  Space  Environment  Service  Center  of  the  National 
Oceanic  and  Atmospheric  Administration  (NOAA)  in  Boulder,  for  the  periods  of  the 
Callaghan  trips. 

Three  solar  phenomena  are  known  to  affect  VLF  propagation: 

1.  Ionization  due  to  X-rays  from  solar  flares  (a  daytime-only  effect). 

2.  Ionization  due  to  Solar  Proton  Events  (SPEs),  which  are  similar  to  Polar-Cap 
Absorption  (PCA)  events  and  mostly  effect  propagation  North  of  60°  N.  latitude. 

3.  Magnetic  field  disturbances  (geomagnetic  storms)  caused  by  the  flow  of  solar 
particles  Into  the  earth’s  magnetosphere,  where  consequent  disturbances  are 
sometimes  severe  enough  to  affect  the  ionosphere.  Severe  geomagnetic  storms 
often  force  the  auroral  zone  South  of  its  normal  60°  N.  position. 

Little  solar  activity  was  anticipated  during  the  1985-1986  period  of  the  Callaghan 
trips,  because  the  solar  cycle  was  near  a  minimum  between  the  usual  eleven-year  peaks 
of  activity.  Only  four  SPEs  occurred  during  the  trip  period,  and  only  three  of  these  co¬ 
incided  with  actual  trip  times.  Table  7  lists  the  start  and  peak  times  (again  in  UT)  for 
applicable  SPEs.  Because  SPE  activity  was  short  in  each  case,  no  correlation  against 
measured  data  effects  was  found. 

Only  five  trips  occurred  during  times  of  significant  geomagnetic  storm  activity.  The 
five  trips,  along  with  relative  levels  of  storm  activity  are  listed  in  table  8.  Storm  levels 
varied  from  "minor”  to  "severe"  in  intensity,  where  NOAA  24-hour  geomagnetic  "A 
Index"  levels  are  defined  as: 

•  Minor  Storm:  A  Index  from  30  to  50 

•  Major  Storm:  A  Index  from  50  to  100 

•  Severe  Storm:  A  Index  >  100 

Only  one  severe  storm  was  experienced  during  the  1985-1986  period,  but  it  occurred 
during  Trip  41,  which  had  been  rejected  from  the  primary  trip  sets  at  all  four  frequencies. 
The  other  storm  levels  were  judged  to  be  too  low  and  infrequent  to  correlate  with  VLF 
measured  data. 
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Table  7.  Solar  proton  events  occurring  during  Callaghan  trips. 


TRIP 

ORIGIN 

DESTINATION 

TRIP  START 

TRIP  END 

SPE  START 

SPE  PEAK 

4 

BAYONNE 

BREMERHAVEN 

04*24-88 

0400 

08*02*88 

1018 

04*28*88 

1430 

04*28*88 

0600 

14 

SAVANNAH 

ALEXANDRIA,  EGYPT 

07*08*88 

0430 

07*19*88 

0230 

07*09*86 

0236 

07*09*88 

0328 

40 

BAYONNE 

CHARLESTON 

02*04*86 

2330 

02*06*88 

1118 

02*08*88 

0828 

02*07*88 

1730 

- 

- 

- 

- 

- 

- 

- 

02*14*66 

1188 

02*18*86 

0400 

Table  8.  Callaghan  trips  during  geomagnetic  storms  (selected  trips  in  large  type). 


STORM 

LEVEL 

TRIP 

ORIGIN 

DESTINATION 

START 

FND 

MAJOR 

4 

BAYONNE 

BREMERHAVEN 

04*24*88 

0400 

08-02*88 

1018 

MINOR 

14 

SAVANNAH 

ALEXANDRIA,  EGYPT 

07-08*88 

0430 

07*19*88 

0230 

MIN./MAJ. 

17 

BAYONNE 

BREMERHAVEN 

08*08*86 

0100 

03*14*88 

2200 

MIN./MAJ. 

37 

JACKSONVILLE,  FL 

ZEEBRUGQE 

12*29*88 

0700 

01*06*86 

1846 

SEVERE 

41 

CHAAuBtTON 

AOHIADAM 

03-04*44 

0400 

03-14-44 

0144 

31 
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II-5.4  PCAP  EVALUATIONS 

Figure  2,  on  page  9,  illustrates  typical  Callaghan  trips  with  respect  to  the  position 
of  the  70°  dip  angle  line  (the  wide  line).  Those  trips  which  pass  North  of  Scotland  are 
seen  to  cross  an  extensive  ocean  area  North  of  the  line.  The  70°  dip  angle  line  is  also 
shown  on  figure  3  with  respect  to  the  four  signal  propagation  paths  selected  for  predicted 
data.  Transmissions  from  the  two  CONUS  sources  at  21.4  and  24.0  kHz  pass  North  of 
the  70°  line  in  the  initial  portion  of  their  paths.  The  transmissions  from  the  British  trans¬ 
mitters,  however,  at  16.0  and  19.0  kHz,  pass  below  the  70®  line. 

The  nighttime  ionosphere  is  known  to  be  more  depressed  in  the  polar  regions,  be¬ 
coming  rrt  re  day-like  at  the  higher  latitudes.  The  LWPC  program  recognizes  this  effect 
by  including  a  polar-cap  (PCAP)  transition  that  modifies  the  /3-h'  model  across  a  user- 
selectable  boundary,  The  default  PCAP  setting  is  70-74,  corresponding  to  a  transition 
area  from  70°  to  74°  in  dip  angle.  /9-h'  profiles  are  automatically  adjusted  for  propa¬ 
gation  paths  passing  through  the  PCAP  transition.  Table  9  is  an  example  listing  of 
segments  along  the  60°  propagation  path  from  the  Annapolis  transmitter,  where  Rho  is 
distance  along  the  patn.  As  the  path  moves  through  the  area  North  of  the  70°  dip  angle 
line,  and  the  codip  angle  falls  below  20°  (codip  -  90®  -  dip),  a  corresponding  change  is 
seen  in  /9-h'.  Table  10  lists  the  Callaghan  trips  that  traversed  areas  below  the  70®  dip 
angle  line,  and  identifies  those  selected  for  the  primary  trip  sets  at  21.4  and  24.0  kHz. 

The  Callaghan  data  is  thus  a  mixed  bag  from  the  standpoint  of  polar-cap  effects,  and 
it  proved  difficult  to  isolate  PCAJP-free  data.  Instead,  the  project  ran  LWPC  in  the  de¬ 
fault  setting  (PCAP -70-74)  for  all  frequencies,  then  added  excursions  at  21.4  and  24.0 
kHz  where  PCAP  was  set  at  80-84,  thus  moving  the  transition  area  well  above  the  prop¬ 
agation  path  areas.  Results  of  the  PCAP  excursions  are  presented  in  the  next  sections. 

Table  9.  Path  segmentation  for  Annapolis,  path  bearing  60®. 

RHO  AZIM  CODIP  MAGFLD  SIGMA  EPSR  BETA  HPRIME 

0.000  068.0  194  0.542  4e  +  00  81.0  0.39  84.8 

0.060  069.1  19,2  0.542  38  -  03  15.0  0.39  84,8 

0.260  072.8  18,7  0,543  4e--00  81.0  0,39  84.8 

0.800  083,1  18.0  0.539  4e-Q0  81.0  0.38  82.7 

1.180  090.0  18.0  0.533  4e-^00  81.0  0.3S  84.B 

2.460  100.5  20.0  0.502  4e-QQ  81.0  0.41  87.0 

5.920  120.0  30.0  0.435  18-02  15.0  0.41  87.0 

6.240  120.8  31.3  0.431  18-02  15.0  0.41  87.0 

6.460  121.4  32.3  0.428  4e-00  81.0  0.41  87.0 

7.060  123.1  35.2  0.421  4e-0Q  81.0  0.41  87.0 
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Table  10.  Callaghan  Trips  below  the  70°  dip  angle  line  (large  type). 


TRIP 

ORIGIN 

DESTINATION 

START 

END 

1 

•BAYQNNI.  1C 

BMMMIAHAVEN 

931119 

HOO 

0401*99 

1700 

2 

•BREMERHAVEN,  W,  OR. 

CHARLESTON 

04.11-85 

2200 

04-19-85 

0900 

3 

CHAAlIITON,  1C 

BAVONNC 

04.19.19 

0900 

04-H-I9 

0700 

4 

-IAYONN1 

9*IM0«HAVIN 

04. 34. IB 

0400 

01-01-19 

1010 

9 

•BHBM1AHAV1N 

BAYONNE 

09-09.I9 

0030 

09-13-19 

0130 

6 

•BAYONNI 

BREMENHAVEN 

OO*. 9*91 

0343 

01-12-01 

0000 

T 

'BREMgftHAVEN 

BAvONNe 

08-33.00 

3100 

09*30*19 

0900 

1 

-1AYONNI 

BRBMIHHAVBN 

00-31  -tO 

8300 

Ot'OMB 

0930 

i 

■MI'.'tfSHAVIN 

BAVONNI 

09-09'IS 

1000 

ot-14-ia 

KOO 

10 

•QAVONNI 

■*|M|MHAV«N 

01-1T-1B 

0400 

01*13*19 

0100 

M 

•IHMiRHAVIN 

BAYONNI 

09-19*99 

0400 

07-01-19 

0730 

■a 

BAVONNi 

CHAStaiTON 

09. 91* IB 

0100 

07-03-16 

0700 

13 

charleston 

SAVANNAH 

07-94.88 

0800 

07-04-88 

0830 

14 

•SAVANNAH 

ALEXANDRIA,  EGYPT 

07.08-88 

0430 

07-18-88 

0230 

IS 

•ALEXANDRIA 

BREMERHAVEN 

07.21-88 

1400 

07-27-86 

2130 

'* 

•»«JM|4HAVIS 

BAYONNI 

071199 

0000 

01-04  19 

1719 

'7 

•BAYONN1 

SABMIflHAVIN 

010919 

0100 

09-14.19 

1100 

IS 

BREMERHAVEN 

ALEXANDRIA 

08-16-88 

0200 

08-22-66 

0330 

19 

•ALEXANDRIA 

SAVANNAH 

08.25*88 

0700 

09-08-85 

0712 

20 

•SAVANNAH,  QA 

BREMERHAVEN 

08-10-88 

0330 

08-18-86 

1900 

ai 

‘1RBMEAHAVIN 

NOMOLK 

09-10-19 

1400 

M-I9-II 

1000 

22 

NORFOLK,  VA 

CHARLESTON 

08-30-85 

0118 

08-30-86 

1800 

23 

NORFOLK 

CHARLESTON 

10-23-88 

1800 

10-24-85 

0600 

24 

•CHARLESTON 

ROTTERDAM,  HOLLAND 

10-28-88 

1300 

11-01-85 

1700 

28 

ROTTERDAM 

LIVERPOOL 

11-03-88 

1448 

11-04-88 

2000 

26 

LIVERPOOL,  ENGLAND 

ANTWERP 

11-08-88 

1834 

11-07-88 

1722 

27 

ANTWERP,  BELGIUM 

BREMERHAVEN 

11-08-66 

2400 

11-08-88 

0830 

28 

•BREMERHAVEN 

CHARLESTON 

11-12-88 

0300 

11-21-85 

0830 

H 

CHAAlllTON 

BAYONNI 

11-11-11 

0100 

11-13-11 

1100 

Jfl 

BAYONNI 

NO"*OIK 

1 1  -11-19 

0110 

11-ll-IB 

1000 

31 

NORFOLK 

SAVANNAH 

12-01-86 

2300 

12-02-68 

1830 

32 

SAVANNAH 

CRISTOBAL  PANAMA 

12-08-85 

1430 

12-08-88 

0841 

33 

CRISTOBAL 

SAVANNAH 

12-10-88 

1740 

12-13-88 

0830 

34 

savannah 

PUERTO  CORTEZ,  HONDURAS 

12-17-88 

0200 

12-19-88 

1122 

39 

NO  OATA 

36 

CRISTOBAL 

JACKSONVILLE 

12-23-85 

2214 

12-28-86 

1100 

37 

•JACKSONVILLE,  FL 

ZSEBRUGC1E 

12-28-88 

0700 

01-06-88 

1848 

38 

Z2E3RUGQE.  BELGIUM 

BREMERHAVEN 

01-00-88 

0854 

oi-oe-ee 

1928 

» 

•bmuirhavin 

BAYONNI 

0M0*M 

1300 

01-lf-M 

1900 

<e 

BAYONNI 

CnAR'JITON 

01 '04-49 

1310 

01-06-91 

1111 

41 

CHARLESTON 

ROTTERDAM 

02-08-86 

0400 

0?.- 16-68 

0184 

49 

NO  DATA 

43 

•BREMERHAVEN 

BEAUMONT 

02-22-88 

0800 

03-04-88 

0718 

44 

BEAUMONT,  TX 

CHARLESTON 

03-07-86 

1830 

03-08-88 

ooia 

45 

•CHARLESTON 

ROTTERDAM 

03-08-86 

0800 

03-16-86 

0600 

<e 

BorinnAM 

IS3IN,  NORWAY 

S3'l-ll 

0300 

03-10-11 

1000 

4’ 

•BOOIN 

MO«|MUDCnY  NC 

0314  91 

0010 

04-01-19 

1300 

48 

MOREHEAD  CITY 

CHARLESTON 

04-08-86 

0018 

04-08-66 

0800 

4® 

ChABlBITON 

IAYONNI 

04-00*10 

0900 

04-07-11 

0130 

■3 

•BAVONNi 

0401-00 

0100 

04-19-11 

MOO 

81 

ROTERDAM 

BREMERHAVEN 

04-18-86 

0230 

04-18-86 

1101 

91 

•B«iMB4HAV*N 

BAVONNI 

04.19-11 

0100 

04-1919 

1700 

•Indicates  trips  selected  lor  Primary  Trip  Sets  at  21.4  and  24,0  kHz. 
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SECTION  II-6 

PREDICTED  VS.  MEASURED  DATA  COMPARISONS 
II-6.1  RUGBY,  16.0  kHz 

Primary  results  of  the  analysis  are  presented  in  this  and  the  following  sections.  As 
was  described  in  the  methodology  discussion  of  section  11*2,  PLOT.DIS  amplitude  vs. 
distance  plots  of  aggregated  data  from  the  Callaghan  recordings  are  compared  with  over¬ 
lays  of  predicted  data  graphed  to  the  same  scale.  Figures  9*38  are  plots  of  such  data 
arranged  as  recorded  and  predicted  pairs,  with  the  measurements  and  predictions  shown 
as  a  function  of  distance  from  the  transmitter  along  great  circle  paths.  All  data  is  for 
nighttime,  extracting  trip  data  satisfying  the  criteria  that  both  the  transmitter  and  re¬ 
ceiver  (ship)  locations  have  solar  zenith  angles  greater  than  90°. 

Graphs  for  the  Rugby  transmitter  at  16.0  kHz  are  given  in  figures  9-13,  which  are 
illustrative  of  the  many  comparisons  made  during  the  analysis.  Because  of  the  great  vol¬ 
ume  of  such  graphs  generated  by  the  project,  only  representative  samples  are  presented 
in  this  report.  A  common  format  is  used  for  all  frequency  sections,  where  the  first  graph 
compares  recorded  data  with  the  LWPC  default-setting  prediction.  A  sequence  of  four 
additional  comparisons  then  follows,  illustrating  a  different  p  selection  spread  over  four 
values  of  h\  The  discussions  of  section  II-7.7  explain  how  summary  graphs  were  used 
to  verify  /S-h'  profile  recommendations  for  the  overall  frequency  range. 

The  16-kHz  examples  of  this  section  indicate  how  the  recorded  vs.  predicted  data 
fit  was  improved  by  increasing  p  and  decreasing  h'  with  respect  to  the  LWPC  default 
settings.  The  comparisons  concentrate  on  the  distance  ranges  between  zero  and  five 
megameters  in  searching  for  best  data  fits.  It  is  possible  to  improve  the  fit  in  the  ex¬ 
treme  far  field  (5-8  Mm)  by  increasing  p  substantially,  but  the  close-in  curves  no  longer 
match  as  well  and  the  selection  becomes  suspect.  Constraint  of  the  data-fit  selection 
to  distances  <5  Mm  is  justified  because,  beyond  that  range,  the  measured  data  is  less 
reliable  as  recorded  levels  approach  those  of  atmospheric  and  shipboard  noise. 

In  each  of  the  figures,  applicable  trip  numbers  are  shown  below  the  lower  edge  of 
the  graph.  The  predicted  data  file  noted  above  the  graph  gives  a  two-character  run  code, 
then  the  letters  “as”  denoting  the  LWPC  ASCII  output  file,  followed  by  the  /3-h'  des¬ 
ignation.  For  example,  the  notation  of  figure  9  shows  that  “Run  fc”  was  for  a  p  range 
of  0.30  to  0.80  (corresponding  to  a  10-  to  60-kHz  frequency  range)  and  was  for  an  h' 
of  87  km. 
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TCS  PLOTJJIS  09/26/89  15:09  Predicted  data  file:  fcas3080.870 

Freq:  16.000  kfreq:  160 

NIGHT  (29-Mar-B5)  to  (25-Apr-86)  deldbl:  8.  deldb2:  12. 


37  39  43  45  47  50  52 

Figure  9.  LWPC  prediction  for  ft = 036,  h'=87  (default),  PCAP= 70-74, 
with  recorded  data  for  Rugby  at  65  kW,  16.0  kHz. 


PLOT _DIS  09/26/89  15c i 3  Predicted  data  file:  beas4030.840 

16.000  kfreq:  160 

HT  (29-Mar -85)  to  (25-Apr-86)  deldbl:  8.  deldb2:  12 


PCAP= 70-74, 
V,  16.0  kHz. 


PLOTJHS  09/26/89  19:  10  Predicted  data  file:  eeas4090.B50 

r  16.000  kfreq:  160 

HI  (29-Mar-85)  to  (25-Apr-86)  deldbl:  9.  deldb2:  12. 


37  39  43  45  47  50  52 

Figure  11.  LWPC  prediction  for  (i — 0.46,  h'=85,  PCAP  =  70-74, 
with  recorded  data  for  Rugby  at  65  kW,  16.0  kHz. 


PLOTJJIS  09/26/89  15:18  Predicted  data  file:  geas4090.B60 

[.  16.000  kfreq:  160 

HT  (29-Mar -85)  to  (25-Apr-86)  deldbl:  8.  deldb2:  12 


70-74, 


PLOTJJIS  09/26/89  15: 22  Predicted  data  file:  feas4Q90.B70 

i:  16.000  kfreq:  160 

HT  (29-Mar-85)  to  (25-Apr-86)  deldbl:  8.  deldb2:  12. 


Di  stance  (Mb) 


I 


11*6.2  ANTHORNE,  19.0  kHz 

Results  for  the  19.0-kHz  transmitter  are  shown  in  figures  14*18.  As  was  the  case 
for  16  kHz,  the  default  profile  shows  the  longer  "cycles"  between  fades  when  compared 
with  measured  data,  indicating  that  a  change  to  a  lower  h'  is  in  order.  This  is  confirmed 
in  figures  15-18,  where  /$  was  adjusted,  also. 

/9  is  increased  to  the  0.40-0.90  range  in  this  instance.  With  interpolation  by  the 
LWPC,  a  /9  of  0.49  was  selected  by  the  computer  for  this  frequency. 


PLOTJHS  09/26/89  15:27  Predicted  data  file:  fcas3080.B70 

|:  19.000  kfreq:  190 

IHT  (29-*4ar-85)  to  (25-Apr-86)  deldbl:  8.  deldb2:  12 


•C  prediction  for  fi  =  039,  h'  =  87  (default),  PCAP= 70-74 
recorded  data  for  Anthonie  at  80  kW,  19.0  kHz. 


TCS  PL0TJ3IS  09/26/09  15:31  Predicted  data  file:  geas4O9O.04O 

Freq:  19.000  kfreq:  190 

•  NIGHT  (29-Mar- 85)  to  (25- Apr -06)  deldbl:  B.  deldb2:  12. 


CD 


a 

?a 


cm 

<n 

CM 

.2 

C  ps. 
(O  <n 
4J 

in  id 

£2  « 

a 

m 


cm 


§$ 

t".  n 

o  v 

8ft 

Sfc 

CM  Ifl 

o  cn 
o 


v 


(ui/An;  aAoqe  gp)  spn^uduiy 
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L  LWPC  prediction  for/? =0.49,  b'=84,  PCAP = 70-74, 
with  recorded  data  for  Anthorne  at  80  kW,  19.0  kHz. 


PL0TJ3IS  09/26/89  15:36  Predicted  data  file:  eeas4090.850 

r  19.0C0  kfreq:  190 

HT  (29-Mar-85)  to  (25-Apr-86)  deldbi:  8.  deldb2:  12 


Anthonie 


PLOTJIIS  09/26/89  15:40  Predicted  data  file:  heas4090.860 

:  19.000  kfreq:  190 

HT  (29-Mar-85)  to  (25-Apr-86)  deldbl:  B.  deldb2:  12 


II-6.3  ANNAPOLIS,  21.4  kHz 

As  frequency  increases,  wavelength  shortens  and  the  earth-ionosphere  waveguide  can 
support  many  more  propagation  modes.  This  is  confirmed  by  the  added  complexity  of 
both  measured  and  predicted  data  at  21.4  kHz.  when  compared  with  the  similar  curves 
of  the  lower  frequencies. 

The  default  LWPC  profile  shown  in  figure  19  reflects  the  addition  of  higher-order 
modes,  but  lacks  the  deep  fading  patterns  suggested  by  the  first  two  nulls  of  the  mea¬ 
sured  data.  These  nulls  are  restored  in  the  £*0.51  profiles  of  figures  20-23,  but  the 
selection  of  the  best  h'  has  become  more  difficult  and  subjective  at  this  higher  frequency. 
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PLOTJJIS  09/26/89  15:50  Predicted  data  file:  icas30B0.870 


C  prediction  tor  fi =0.41,  h'=87  (default),  PCAP  =  70-74. 
recorded  data  for  Annapolis  at  250  kW,  21.4  kHz. 


PLOT J3IS  09/26/89  16:03  Predicted  data  file:  geas4090.860 

F  21.400  kfreq:  214 

HT  (30-Mar -85)  to  (25-Apr-86)  deldbl:  8.  deldb2:  12 


“C  prediction  for  /J=0-51,  h'=86,  PC AP— 70-74. 
recorded  data  for  Annapolis  at  250  kW,  21.4  kl 


PL0TJ3IS  09/26/89  19:23  Predicted  data  file:  ieas4090.B70 

i:  21.400  kfreq:  214 

MT  (30-Mar -85]  to  (25-Apr-B6)  deldbl:  8.  deldb2:  12 


LWPC  prediction  for  £=0.51,  h'=87,  PCAP= 70-74, 
with  recorded  data  for  Annapolis  at  250  kW,  21.4  kl 


11*6.4  CUTLER,  24.0  kHz 

The  default  LWPC  prediction  of  figure  24  shows  even  more  high-order  mode  effects 
than  in  the  21.4-kHz  case.  Again,  the  nulls  seem  less  pronounced  than  in  the  measured 
data. 

Selection  of  a  higher  fi  range,  as  in  figures  25*28,  results  in  a  partial  improvement, 
but  the  match  of  fading  positions  is  disappointing.  This  sequence  suggests  the  need  for 
automation  of  the  selection  process  so  that  many  more  profiles  may  be  examined  within 
practical  time  constraints.  Mechanization  of  the  comparisons  is  discussed  in  section  II* 
7.7, 

The  predicted  curves  shown  thus  far  in  the  figures  utilized  the  LWPC  default  polar- 
cap  settings  (see  discussion  in  section  11*5.4).  As  a  further  point  of  interest,  PCAP 
adjustments  were  eliminated  in  the  excursions  of  the  following  sections. 
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CS  PLOTJJiC  09/26/B9  19:55  Predicted  data  file:  keas409C.B40 

ireq:  24.000  kfreq:  240 

NIGHT  (29-J4ar-05)  to  (25-Apr-BB)  deldbl:  8.  deldb£  12. 


Oistam 


CS  PL0TJ3IS  09/26/89  20:00  Predicted  data  file:  feas4Q90.B50 

req:  24.000  kfreq:  240 

NIGHT  (29-Mar-85)  to  (25-Apr-86)  deldbl:  0.  deldb2:  12. 


»  21  24  28  37 


CS  PLOTJJIS  09/26/89  20:09  Predicted  data  file:  ieas4090.870 

ireq:  24.000  kfreq:  240 

NIGHT  (29-Har-85)  to  (25-Apr-86)  deldbl:  8.  deldb2  12. 


PCAP= 70-74, 
kW,  24.0  kHz. 


11*6.5  EXCURSIONS,  ANNAPOLIS,  21.4  kHz,  WITH  PCAP -80*84 

The  graphs  of  figures  29-33  repeat  the  profile  selections  of  section  II-6.3,  except  for 
the  elimination  of  the  LWPC’s  polar-area  adjustment  of /3-h'  by  setting  PCAP -80*84. 
Comparison  of  the  two  sets  of  graphs  shows  the  introduction  of  even  more  high-order 
modes,  and  the  best-fit  selection  is  still  a  difficult  choice. 

One  trend  is  evident:  For  this  selection  of  0,  deletion  of  the  PCAP  adjustment  per¬ 
mits  a  better  match  to  measured  data  at  a  lower  h\  This  group  of  graphs  favors  the 
curve  fit  of  h'-85,  whereas  in  the  PCAP  -70*74  cases  a  higher  h'  is  preferred. 
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PL0TJ3IS  09/26/89  19:28  Predicted  data  file:  deas3080.870 

i:  21.400  kfreq:  214 

;HT  (30-Mar-85)  to  (25-Apr-86)  deldbi:  8.  deldb2:  12. 
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39  43  45  47  50  52 

Figure  29.  LWPC  prediction  for  /?=0.41,  h'=87,  (default),  PCAP = 80-84, 
with  recorded  data  for  Annapolis  at  250  kW,  21.4  kHz. 


PLOT J)IS  09/26/89  19:32  Predicted  data  file:  laas4090.840 

t:  21.400  kfreq:  214 

HT  (30-Mar-85)  to  (25-Apr-88)  deldbl:  8.  deldb2:  12. 


Figure  30.  LWPC  prediction  for  £  =  0.51,  h'  =  84,  PCAP= 80-84, 

with  recorded  data  for  Annapolis  at  250  kW,  21.4  kHz. 


CS  PL0TJ3IS  09/26/89  19:37  Predicted  data  file:  eeas4090.850 

req:  21.400  kfreq:  214 

NIGHT  (30 -Mai — 85)  to  (25-Apr-B6)  deldbl:  6.  deldb2:  12. 


sijodeunv 


CS  PLOTJHS  09/26/89  19:41  Predicted  data  file:  lbas4090.B60 

req:  21.400  kfreq:  214 

NIGHT  (30-Mar -85)  to  (25-Apr-86)  deldbl:  8.  deldb2:  12. 


Distance  (Ha) 


PLOT _pis  09/26/89  19:46  Predicted  data  file:  lcas4090.870 

21.400  kfreq:  214 

MT  (30-Mar-85)  to  (25-Apr-86i  deldbl:  B  deldb2:  12 


Distance  (Hb) 


11-6,6  EXCURSIONS,  CUTLER,  24.0  kHz,  WITH  PCAP-  80-84 

Similar  tradeoffs  were  noted  in  the  24.0-kHz  cases  of  figures  34-38.  The  selection 
of  a  higher  /?  range  and  the  trend  to  lower  h'  values  resulted  in  better  curve  fits  than 
in  the  default  PCAP  graphs. 

No  final  conclusions  could  be  reached  about  the  PCAP  adjustment  in  these  limited 
tests;  data  fit  seemed  to  suggest  a  tradeoff  of  PCAP  range  vs.  h'.  NOSC  has  suggested 
that  more  extensive  tests  of  PCAP  settings,  taken  over  a  much  wider  ^-h'  range,  should 
be  conducted,  and  have  recommended  at  least  three  ranges  for  the  PCAP  transition  anal¬ 
ysis:  80-84,  70-74  and  66-70. 
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PL0TJ3IS  09/26/89  21:35  Predicted  data  file:  laas4090.840 

:  24.000  kfreq:  240 

HT  (29-Mar- 85)  to  (25-Apr-8B)  deldbl:  9.  deldb2:  12. 


CS  PLOT  JUS  09/26/89  21:40  Predicted  data  file:  eeas4090.850 

Ireq:  24.000  kfreq:  240 

WIGHT  (29- Mar-85)  to  (25-Apr-06)  deldbl:  B.  deldb2:  12. 


Figure  36.  LWPC  prediction  for  fi — 034,  b  =85,  PCAP=80-84, 
with  recorded  data  for  Cutler  at  1000  kW,  24.0  kHz. 


CS  PLOT _DIS  09/26/89  21:49  Predicted  data  file:  lcas4090.870 

req:  24.000  kfreq:  240 

NIGHT  (29-Mar-85)  to  (25-Apr-86)  deldbl:  8  deldb2:  12. 


Distam 


SECTION  II-7 
STATISTICAL  ANALYSES 


II-7.1  RUGBY,  16.0  kHz 

While  the  graphs  of  the  previous  section  ore  useful  for  quick  comparisons  of  the  total 
data  set  by  the  human  eye,  the  NOSC  program  TR1PSTAT  provides  a  better  technique 
for  quantitative  evaluation.  In  section  11*7,  statistical  graphs  from  TRIPSTAT  runs  are 
presented  for  each  case  of  the.  earlier  sequence. 

The  statist^al  charts  of  recorded-minus-predicted  signal  strength  in  dB  at  16.0  kHz 
are  given  in  figures  39-43.  Each  chart  presents  separate  histograms  at  125-km  intervals 
along  the  path,  evaluating  the  distribution  of  data  points  recorded  within  a  range  of  ±60 
km  from  each  125-km  interval.  The  total  number  of  recorded  points  within  each  ±60* 
km  range  are  graphed  in  the  lower  chart  of  each  sheet. 

Each  histogram  is  divided  into  1-dB  bins,  the  height  of  which  varies  in  proportion 
to  the  number  of  data  points  in  each  1-dB  interval.  In  addition,  mean  (diamond)  and 
90%  exceedance  level  (x)  values  are  displayed  below  each  histogram.  An  overall  mea¬ 
sure  of  goodness  is  the  degree  with  which  the  histogram  and  its  mean  are  centered  on 
the  0-dB  axis.  A  0-dB  mean  with  all  data  points  centered  on  the  axis  would  represent 
a  perfect  match  between  recorded  and  predicted  signal  levels  at  that  distance. 

The  statistical  graph  for  the  LWPC  default  profile  at  16  kHz  is  shown  in  figure  39, 
where  the  cyclic  nature  of  the  recorded-minus-predicted  data  reflects  the  mismatch  be¬ 
tween  the  fading  cycles  displayed  in  figure  9,  page  35.  The  distribution  of  means  is 
Improved  somewhat  by  the  series  of  charts  in  figures  40-43,  but  these  quantitative  anal¬ 
yses  confirm  that  the  nighttime  ionosphere  is  quite  variable  with  respect  to  the  much 
closer  fits  of  daytime  data  in  the  earlier  Callaghan  report.  This  quantitative  evaluation 
nevertheless  permits  a  better  evaluation  of  one  profile  against  the  other,  as  is  seen  by 
simply  scanning  the  distribution  of  mean  values. 

The  statistical  graphs  were  the  main  means  for  optimizing  the  selection  of  best-fit 
profiles  in  this  task,  considering  the  limited  number  of  profiles  examined.  A  method 
for  automating  this  selection  over  a  greater  number  of  choices  is  discussed  in  section 
II-7.7. 
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TCS  TRIPSTAT  09/18/89  16:  10 

Predicted  data  file:  fcas3080.870 

Freq:  16.000  Months:  01  to  12  MIGHT  Deldis:  60  Xmax: 


Figure  39.  Recorded  -  predicted  data,  Rugby  at  65  kW,  16.0  kHz, 

with  default  LWPC  prediction  for  fi=036,  h'  =  87,  PCAP  =  70-74. 


II-7.2  ANTHORNE,  19.0  kHz 

The  statistical  comparison  between  the  LWPC  default  and  the  optimized  /9-h'  values 
is  more  favorable  in  the  19.0*kHz  examples  of  figures  44-48.  The  distribution  of  means 
in  the  default  case  is  scattered  over  a  wide  range,  while  the  means  of  the  85-km  example 
of  figure  46  are  relatively  well  behaved,  with  most  falling  within  a  ±4-dB  range  out  to 
4.5-Mm  distance. 

The  distribution  of  90%  exceedance  values  is  seen  as  a  less  favorable  indicator  of 
fit  than  is  that  of  the  means  in  figure  46,  where  the  existence  of  spurious  data  points, 
particularly  around  fades  as  is  typical  of  the  nighttime  ionosphere,  produces  large  values 
for  the  90%  exceedance  points. 
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TCS  TRIPSTAT  09/17/89  17:  38 

Predicted  data  file:  geas4090.B40 

Freq:  19.000  Months:  01  to  12  NIGHT  Deldis: 


Figure  45.  Recorded  —  predicted  data,  Anthonie  at  80  kW,  19.0  kHz, 
with  LWPC  prediction  for/?  =  0.49,  h'=84,  PCAP  =  70-74. 


TCS  THIPSTAT  09/21/89  20:50 

Predicted  data  file:  heas4090.880 

Freq:  19.000  Months:  01  to  12  NIGHT  Deldis:  60  Xmax: 


Figure  47.  Recorded  —  predicted  data,  Anthorne  at  80  kW,  19.0  kHz 
with  LWPC  prediction  for  5=0.49,  b*  =86,  PCAP= 70-74. 


48.  Recorded  -  predicted  data,  Anthonie 
with  LWPC  prediction  for  ft — 0.49,  h': 


II-7.3  ANNAPOLIS,  21.4  kHz 

The  LWPC  default  profile  provides  a  much  better  fit  in  the  21.4-kHz  case,  as  is  seen 
by  the  charts  of  figure  49.  The  charts  of  figures  50*53,  for  the  “0.40*0.90  profile  are 
worse  than  those  of  figure  49  at  the  shorter  waveguide  heights.  The  best  fit  for  this 
profile  sequence  is  found  for  h'  “87  km,  in  figure  53,  where  all  but  one  mean  value  fall 
within  a  ±3-dB  range,  and  this  fit  extends  past  the  5-Mm  distance. 
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(gp)  pa^pajd  -  papjooaa  a^unoo 


87 


Distance  (Mb) 

rips:  01  02  04  05  06  07  08  09  10  11  14  15  16  17  19  20  21  24  28  37 
39  43  45  47  50  52 

Figure  53.  Recorded  -  predicted  data,  Annapolis  at  250  kW,  21.4  kHz, 
with  LWPC  prediction  for  /J=0.51,  h'=S7,  PCAP  =  70-74. 


II-7.4  CUT1  SR,  24.0  kHz 

The  statistical  evaluations  for  24.0  kHz  show  that  the  /9-h'  -0.40-0.90  profile  offers 
an  improvement  over  the  LWPC  default  profile,  as  is  shown  in  figures  54*58.  The  dis¬ 
tribution  of  mean  values  is  again  best  at  h'»87  km,  but  fewer  points  are  within  the 
±3-dB  range. 
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TCS  THIPSTAT  09/1B/B9 

Predicted  data  file:  icas30B0.B70 
Freq:  24.000  Months:  01  to 
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Figure  54.  Recorded  —  predicted  data.  Cutler  at  1000  kW,  24.0  kHz, 

with  default  LWPC  prediction  for  0=0.44,  h'  =  87,  PCAP  =  70-74. 


TCS  TRIPSTAT  09/17/89  19:29 

Predicted  data  file:  keas4090.B40 

Freq:  24.000  Months:  01  to  12  NIGHT  Deldis:  60  x«ax: 


Recorded  -  predicted  data.  Cutler  at  1000  kW,  24.0  kHz, 
with  LWPC  prediction  for  fi =0.54,  h'  =  84,  PCAP  =  70-74. 


TCS  THIPSTAT  09/17/09  19:00 

Predicted  data  file:  ieas4090.870 


Recorded  —  predicted  data.  Cutler  at  1000  kW,  24.0  kHz, 
onth  I  WPT  nmiictinn  for  B= 054.  b'  =87.  PCAP  =70-74. 


II-7.5  EXCURSIONS,  ANNAPOLIS,  21.4  kHz,  WITH  PCAP -80-84 

Figures  59-63  examine  the  statistical  characteristics  of  the  21.4-kHz  path  with  the 
PCAP  adjustment  of  LWPC  disabled.  As  was  seen  in  the  signal  vs.  distance  plots,  re¬ 
moval  of  the  PCAP  compensation  forces  the  best-fit  case  to  a  lower  reflection  height. 
The  best  case  for  these  examples  is  found  at  h'-85  km. 

In  the  85-km  charts  of  figure  61,  most  mean  values  to  beyond  5-Mm  distance  are 
within  ±2  dB  of  the  zero  line.  This  agreement  between  measured  and  predicted  data 
is  indeed  remarkable,  given  the  variable  nature  of  the  nighttime  ionosphere,  and  ap¬ 
proaches  the  agreement  achieved  for  the  daytime  cases  of  the  Callaghan  report. 
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TCS  TPIPSTAT  09/21/09  22:  02 

Predicted  data  file:  lbas4090.B60 

Frea:  21.400  Months:  01  to  12  NIGHT  Deldis 
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Figure  62.  Recorded  -  predicted  data,  Annapolis  at  250  kW,  21.4  kHz, 
with  LWPC  prediction  for  /J=l).5l,  h'  =86,  PCAP=80-84. 


II-7.6  EXCURSIONS,  CUTLER,  24.0  kHz,  WITH  PCAP*  80-84 

The  excursions  of  figures  64-68,  for  24.0  kHz,  reflect  the  variability  of  the  signal  vs. 
distance  plots  seen  earlier.  Although  the  fit  isn’t  as  good  as  in  the  21.4-kHz  examples, 
the  best  case  is  still  found  at  84  to  85  km.  This  suggests  that  the  85-km  profile  be  se¬ 
lectee  as  a  standard  improvement  for  all  nighttime  frequencies,  since  it  would  then  be 
a  common  choice.  (A  preliminary  recommendation  for  nighttime  profiles  from  the  com¬ 
parisons  of  this  task  was  presented  in  figure  1  on  page  5.)  However,  such  a  decision 
ignores  the  effect  of  LWPC’s  PCAP  adjustment,  which  certainly  must  remain  in  place 
at  some  setting  to  compensate  for  the  lower  ionospheres  at  higher  latitudes. 

These  findings  underscore  the  need  for  additional  analyses  of  the  Callaghan  data  over 
a  wider  range  of  profiles,  and  for  the  investigation  of  more  wide-ranging  PCAP  selec¬ 
tions.  The  manner  in  which  the  statistical  comparisons  could  be  automated  is  discussed 
in  next  section. 
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TCS  THIPSTAT  09/18/89  17:34 

Predicted  data  file:  deas30B0.870 

Free-  24.000  Months:  01  to  12  NIGHT  Deldis:  60  Xmax 


Figure  64.  Recorded  —  predicted  data.  Cutler  at  1000  kW,  24.0  kHz, 

«  h  default  LWPC  prediction  for  /?=0.44,  h'  =87,  PCAP  =  80-84. 


TCS  TRIPSTAT  09/22/89  09:00 

Predicted  data  file:  lbas4090.860 

Freq:  24.000  Months:  01  to  12  NIGHT  Deldis: 


m 

8 

««r 

ru 

•n 

ru 

S3 

|S 

w  K 

0)  «H 

U 

c  in 
(tJ  Jh 

to  in 


s 


ru 


r»  oj 
o  in 

U0  Q 

o  in 
in 

o  n 

m 

o 

oj  m 


«  cn 
o  n 


to 

a 


wnon 


Figure  67.  Recorded  —  predicted  data.  Cutler  at  1000  kW,  24.0  kHz, 
with  LWPC  prediction  for  0=0.54,  h'  =86,  PCAP=80-84. 


TCS  THIPSTAT  09/22/89  09;  18 

Predicted  data  file:  lcas4090.870 

Freq;  24.000  Months:  01  to  12  MIGHT  Deldis: 


(SP)  pe*3TP8Jd  -  papjooey  e^unoo 
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Recorded  —  predicted  data.  Cutler  at  1000  kW,  24.0  kHz, 
with  LWPC  prediction  for  /J=(t54,  b'=87,  PCAP=80-84. 


II-7.7  SUMMARY  OF  STATISTICAL  COMPARISONS 

The  previous  sections  showed  that  TRIPSTAT  statistical  charts  provide  a  means  for 
quantitative  evaluation  of  recorded-minus-predicted  data  that  is  superior  to  the  analyst's 
subjective  judgment,  Still,  the  comparison  of  one  statistical  chart  against  another  re¬ 
quires  a  similar  visual  appraisal. 

As  a  way  of  further  mechanizing  the  comparisons,  an  average  of  all  the  histogram 
mean  values  was  determined  for  each  chart  over  given  distance  ranges.  Then,  a  single 
number  was  derived  for  each  profile  and  distance  range,  and  these  numbers  were 
graphed  over  a  range  of  profile  choices.  Graphs  of  the  average  of  mean  values  appear 
in  this  section,  and  represent  summaries  of  the  final  profiles  selected  for  comparison. 

Similarly,  standard  deviations  of  each  histogram  were  computed  from  the  mean  and 
90%  exceedance  levels  and  were  averaged  for  each  chart  over  a  particular  distance 
range,  assuming  a  normal  distribution  of  data.  Again,  a  single  number  corresponding 
to  average  standard  deviation  was  derived  for  each  profile  and  distance,  and  these  num¬ 
bers  were  graphed  for  a  range  of  profile  choices. 

Figures  69-74  present  graphs  of  averages  of  the  mean  and  standard  deviation  values 
vs.  h'  for  the  selected  0  range  of  0.40-0.90.  Figures  69-72  summarize  the  four  analysis 
frequencies,  using  the  default  PCAP  settings  of  70-74,  Summaries  of  the  21.4-  and  24.0* 
kHz  excursions  with  PCAP  ■  80-84  are  given  in  figures  73-74.  In  each  graph,  averages 
for  five  distance  ranges  are  given,  and  a  point  representing  the  single  average  for  the 
LWPC  default  profile  at  a  0-5  Mm  range  is  also  shown. 

Choosing  a  profile  where  the  summary  graph  average  of  means  is  at  a  minimum  is 
the  final  step  of  the  analysis  process.  For  example,  in  figure  70  for  19-kHz,  the  0-5  Mm 
curve  shows  the  lowest  average  mean  at  an  h'  just  below  85  km. 

A  scan  of  the  summary  curves  confirms  that  h'»85  km  is  a  good  choice  for  the 
(3  **  0.40-0.90  profile  range,  providing  that  PCAP  ■  80-84  at  the  higher  frequencies.  Ide¬ 
ally,  three-dimensional  plots  with  both  /3  and  h'  treated  parametrically  over  a  much  wider 
profile  range  should  be  produced,  and  PCAP  adjustment  effects  should  be  investigated. 

The  summary  graphs  show  how  the  technique  for  finding  the  average  of  means  and 
standard  deviations  can  aid  in  the  optimization  process.  The  summaries  of  this  task  were 
derived  manually  from  the  TRIPSTAT  charts.  Future  tasks  should  facilitate  the  process 
by  adding  program  routines  to  automatically  produce  listings  and  graphs  of  such  averages 
from  a  full  range  of  0-h'  profiles. 
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Average  of  Sid.  Dev.  Values  (dB)  Average  of  Ueon  Values  (dB) 


83.5  84  84.3  85  85.5  86  86.5  87  87.3 

h  *  (km) 

Dtitance 

Range  (Mm):  •  0-1  ••  0-2  —  0-3  —  0-4  —  0-5 

e  For  Beta  0.3-0.8,  h'  87  (LWPC  Defoult)  at  0-5  Mm. 


Figure  69.  Average  of  mean  and  standard  deviation  values  vs.  h\ 
recorded  -  predicted  data  for  Rugby,  65  kW,  16.0  kHz, 
for  p  range  0.40-0.90,  PCAP- 70-74. 
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tkiswoxht  Average  of  Mean  Values 


83. S  84  84.5  85  85.5  86  86.5  87  67.5 

h  '  (km) 

Average  of  Standard  Deviation  Values 


63.5  84  84.5  85  85.5  86  86.5  87  87.5 


h  *  (km) 

Distance 

Rang#  (Mm):  0-1  -•  0-2  -0-3  —  0-4  —0-5 

9  For  Bata  0. 3-0.8,  h*  87  (LWPC  Dafault)  at  0-5  Mm. 


Figure  70.  Average  of  mean  and  standard  deviation  values  vs.  h', 

recorded  -  predicted  data  for  Anthorne,  80  kW,  19.0  kHz, 
for  range  0.40-0.90,  PCAP  ■  70-74. 
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Average  of  Std.  Dev.  Values  (dB)  Average  of  Mean  Values  (dB) 


1 


h  *  (km) 

Distance 

Rang#  (Mm):  0-1  -■  0-2  —  0-3  —  0-4  — —  0-5 


e  For  B«ta  0.3-0.8,  h'  87  (LWPC  Default)  at  0-5  Mm. 


Figure  72.  Average  of  mean  and  standard  deviation  values  vs.  h', 

recorded  -  predicted  data  for  Cutler,  1000  kW,  24.0  kHz, 
for  /9  range  0.40-0.90,  PCAP«  70-74. 


110 


Average  of  Std.  Dev.  Values  (dB)  Average  of  Mean  Values  (dB) 
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Average  of  Standard  Davlatlon  Values 
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0 4 — ‘  -i  — i —  >■—  i — ■  "" i  ■  *■  ■■■■■*  ■■■  ■  t  1  I 

83.5  84  84.5  85  85.5  86  66.5  B7  87.5 

h  •  (km) 

Distance 

Range  (Mm):  •  0-1  --  0-2  —  0-3  —  0-4  — —  0-5 

©  For  Beta  0.3-0.8,  h*  87  (LWPC  Default)  at  0-5  Mm. 


Figure  73.  Average  of  mean  and  standard  deviation  values  vs.  h\ 

recorded  -  predicted  data  for  Annapolis,  250  kW,  21.4  kHz, 
for  range  0.40-0.90,  PCAP*  80-84. 


Ill 


Average  of  Sid.  Dev.  Values  (dB)  Average  of  Mean  Values  (dB) 


h  *  (km) 

Distance 

Range  (Mm):  0-1  -•  0-2  -0-3  —  0-4  —  0-5 

q.  For  Beta  0.3-0. 8,  h*  87  (LWPC  Default)  at  0-5  Mm. 


Figure  74.  Average  of  mean  and  standard  deviation  values  vs.  h', 

recorded  -  predicted  data  for  Cutler,  1000  kW,  24.0  kHz, 
for£  range  0.40-0.90,  PCAP*  80-84. 
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SECTION  8 


CONCLUSIONS 

All  objectives  of  Part  1  of  the  present  task  have  been  met,  and  the  following  con¬ 
clusions  have  been  made. 

•  The  analysis  methodology  of  the  Callaghan  report  has  been  repro¬ 
duced,  and  comparable  output  graphics  and  data  reduction 
techniques  have  been  reduced  to  practice.  The  NOSC  LongWave 
Prediction  Capability  (LWPC)  family  of  programs  has  been  installed 
on  task  computers  in  an  emulated  VAX- VMS  environment  and  us¬ 
age  has  been  validated  against  NOSC-supplied  test  cases. 

•  The  daytime  analyses  of  the  Callaghan  report  have  been  repeated 
with  nearly  identical  results,  as  a  further  measure  of  methodology 
validation.  Example  predictions  of  the  same  type  have  been  pro¬ 
duced  for  the  nighttime  environment,  and  initial  comparisons  with 
recorded  data  at  21.4  kHz  show  good  agreement. 

•  Statistical  summary  graphs  should  be  expanded  to  also  include  dis¬ 
tance  ranges  of  1-5,  2-5,  3-5  and  4-5  Mm. 

•  It  is  evident  that  separate  graphs  of  recorded  data  by  trip  will  be 
helpful  in  eliminating  questionable  data  and  in  optimizing  a  pre¬ 
ferred  set  of  aggregate  data  for  each  frequency.  Such  optimized  sets 
of  recorded  data  should  be  routinely  assembled  before  profile  mod¬ 
ifications  are  attempted. 

•  Because  the  presence  of  significant  high-order  mode  interaction 
tends  to  modify  the  dominant  mode  propagation  signature,  and  be¬ 
cause  this  effect  varies  with  frequency,  the  optimization  of  nighttime 
profiles  is  expected  to  be  complex.  Profile  modifications  for  the 
nighttime  VLF  environment -the  main  objective  of  this  task- 
should  therefore  be  recommended  only  after  a,l  available 
frequencies  are  considered. 

•  The  Callaghan  report  examined  four  VLF  frequencies:  16.0,  19.0, 

21.4  and  24.0  kHz.  Data  from  the  other  VLF  channels  of  the  re¬ 
corded  database  (16.4,  24.8  and  28.5  kHz)  should  also  be  examined 
within  the  constraints  of  time  and  funding. 

•  If  possible,  future  analyses  should  separate  recorded  data  into  areas 
of  latitude  below  and  above  the  codip  =  20°  point  where  the  present 
LWPC  design  forces  a  profile  change. 
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